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1.0  INTRODUCTION 


During  tlie  past  five  years  an  ONR-sponsored  research  program  has 
been  underway  in  ttie  Klectrical  Kngineering  Department  at  North  Carolina 
State  University  whicli  basically  involves  the  use  of  the  Monte  Carlo 
metliod  to  study  hlgh-fleld  traiisp  .  t  properties  of  semiconductor  nuiterials. 
The  emphasis  has  been  directed  toward  the  study  of  ternary  and  quaternary 


III-V  compounds  in  an  attempt  to  identify  specific  materials  with  desirable 
electronic  properties  for  microwave  and  high-speed  device  applications. 

A  substantial  effort  has  been  devoted  to  developing  a  comprehensive 
Monte  Carlo  computer  program  which  can  be  used  as  a  predictive  tool  in 
the  study  of  hot  electron  transport,  as  well  as  low  field  ohmic  transport, 
in  a  general  class  of  semiconductor  materials.  We  presently  believe  that 
the  capabilities  which  have  been  developed  so  far  under  ONR  sponsorship 
are  second  to  none  within  the  scientific  community,  and  that  the  bulk 
static  transport  properties  of  two  general  classes  of  quaternary  III-V 
semiconductors  can  be  predicted  using  our  computational  methods.  These 


two  classes  of  quaternary  materials  consist  of  compounds  of  the  form 

III„III  V  V  III„II1  II1,V  ,  ,  V  V, 

A,  B  C,  D  and  A  B  C  D  (or  conversely  A  BCD),  which  can 
1-xxl-yy  xyz  xyz 

be  specified  to  include  any  binary  or  ternary  compound  as  a  subset  of 


these  general  material  forms. 


The  accuracy  of  the  computational  program  depends  on  several  factors. 


One  major  factor  is  the  way  in  which  the  quaternary  material  parameters 
are  calculated.  The  technique  depends  on  a  knowledge  of  the  binary 
material  parameters  and  either  theoretical  models  or  experimental  data 


for  the  ternary  material  parameters.  This  information  is  coupled  into 


a  general  interpolation  procedure  for  calculating  tlie  quaternary  material 


parameters,  from  wliich  the  transport  properties  are  obtained  in  the 


Monte  Carlo  metliod.  This  work  resulted  over 


prediction  that  specific  compositions  of  the  quaternary  compound 


Ga._  In  1’  As  lattice-matched  to  Inl’  substrates  possess  transport 


properties  having  distinct  advantages  over  GaAs  [1]  for  device  applications 


In  the  past  three  years  considerable  experimental  work  on  this  compound  has 


shown  the  material  to  indeed  be  promising  for  many  device  applications 


although  the  principal  applications  have  been  for  electrooptic  devices 


However,  the  complete  advantages  due  to  the  predicted  superior  trans 


port  properties  have  not  been  realized  at  the  present  time,  and  efforts 


to  make  a  microwave  MESFET  have  only  been  partially  successful  [2],  It 


appears  that  many  of  the  device  problems  can  be  overcome  through  the 


refinement  of  the  materials  technology,  and  there  are  reasons  to  remain 


optimistic  that  the  predicted  material  advantages  will  be  realized 


In  the  authors’  opinion  the  influence  of  alloy  scattering  still  remains 


an  unanswered  question.  However,  recent  results  on  Ga  A1  As  and 


Ga  In  As  suggest  that  our  approach  which  uses  the  electron  affinity 


to  be  noted  that  these  recent  results  have  come  about  through  careful 


refinement  of  the  growth  process  in  these  materials  to  the  point  where 


reasonably  pure  material  can  be  achieved.  Then,  by  examining  the  low 


temperature  (77K)  Hall  mobility,  the  effects  of  alloy  scattering  in 


Ga.  A1  As  can  be  separated  from  those  due  to  ionized  impurity  scattering 


and  are  in  good  agreement  with  our  previous  models 
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In  rnntrasL  In  t  linso  rosiilts  vnry  I  i  1 1 1  n  l).isli'  1  iit  ii  rniat  1  on  lias  been 

.ic  cumulated  wliicli  elucidates  exjuT  iinenta]  1  y  the  eli-ctron  transport 

properties  ot  Oa  In  1’  As  which  will  direcllv  determine  its  uselulness 
1  -X  X  1 -y  y 

as  a  material  tor  high-1  reijuency,  high-sp'ed  device  app  1  icat  Ions.  I'rom 

our  own  laboratory  work  in  tlie  growtli  of  the  ciuatorn.irv  it  appears  that 

iiuiny  material  problems  remain  to  be  solved  before  tlie  true  nature  ol  the 

transport  properties  can  be  elucidated.  I'or  example,  energy  bandgap 

measurements  do  support  the  interpolation  procedure  for  determining 

quaternary  material  energy  bandgaps  [A-6J.  However,  measurements  of 

effective  mass  by  various  tecliniques  are  contradictory.  Some  of  the 

measurements  support  the  interpolation  [6,7]  while  other  experimental 

data  are  contradictory  to  it  [8,9].  In  addition,  Hall  effect  measurements 

of  the  temperature  dependence  of  the  electron  mobility  [10,11]  indicate 

that  alloy  scattering  is  considerably  more  prominent  than  is  predicted  by  even 

the  largest  of  the  theoretical  estimates.  In  addition,  for  the  quaternary 

compositions  lattice-matched  to  InP  it  has  been  very  difficult  to  form 

Schottky  barriers  on  n-  or  p-type  material  [2],  and  there  are  often 

problems  in  forming  rectifying  p-n  junctions  by  LPE  [12].  The  emerging 

data  strongly  suggest  regions  of  non-stoichlometric  or  conducting-type 

defects  which  affect  the  barrier  region.  Recent  experiments  in  InP 

grown  by  MBE  [13, lA]  show  the  presence  of  conducting  In  precipitates 

which  affect  the  transport  properties  and  the  formation  of  Au-Schottky 

barriers  on  InP  in  a  manner  which  is  similar  to  the  measurements  on 

Ga,  In  P,  As  [15].  The  question  of  the  behavior  of  impurities,  the 
1-x  X  1-y  y 

formation  of  impurity-defect  complexes,  and  the  behavior  of  non-stoichlo- 
mctrlc  defects  is  emerging  as  a  significant  research  problem  in  the 


III-V  materials  which  will  greatly  affect  the  experimental  and  theoretical 
interpretat ion  of  the  transport  properties.  We  are  beginning  to  examine 
ways  to  study  such  problems  with  the  Monte  Carlo  method,  for  example  by 
initially  including  the  Seitz  tlisloeatlon  scattering  in  the  calculation  [16]. 
The  general  purpose  of  tlie  .ibove  discussion  is  to  show  that  at  the  present 
time  there  is  not  adequate  evidence  to  either  confirm  or  disprove  the 
basic  interpolation  procedure  used  for  predicting  ternary  and  quaternary 
material  parameters.  In  fact,  as  experimental  information  continues  to 
be  gathered,  the  interpolation  procedure  can  be  refined  to  make  it  more 
generally  applicable. 

We  have  continued  during  this  year  to  examine  the  general  area  of 
submicron  device  physics,  and  have  emphasized  the  effort  to  establish 
physical  processes  which  effect  the  operation  of  small  devices  and  to 
incorporate  such  mechanisms  into  the  Monte  Carlo  method.  We  are  examining 
the  limits  of  the  Monte  Carlo  method,  and  other  transport  techniques,  as 
to  their  general  applicability  to  small  device  phenomena. 

As  a  first  step  in  this  approach  we  have  implemented  a  conventional 
transient  or  step-field  approach  [17]  to  examining  what  is  commonly  called 
velocity  overshoot  in  the  electron  dynamics  (18].  Another  approach, 
called  the  "ensemble"  Monte  Carlo  method,  has  also  been  implemented  which 
critically  examines  various  estimation  schemes  for  the  dynamical  behavior 
of  both  electron  drift  velocity  and  electron  diffusion.  We  are  continuing 
to  consider  and  investigate  other  areas  such  as  fluctuation  phenomena, 
modifications  of  collision  interactions  and  interaction  time  field  effects, 
two-dimensional  material  and  device  effects,  fundamental  frequency  limita¬ 
tions,  and  the  Influence  of  degenerate  carrier  statistics. 


Ill  addition,  wo  iiave  bej;iin  to  cxainino  electron  transport  in  quasi- 
two  dimensional  systems  or  supi'r  1  at  t  i  ee  structures.  We  believe  tliat  tills 
emerging  area  has  important  applications  lor  new  device  (.oncepts  and 
structures,  and  that  the  basic  electron  tr  insport  properties  can  be 
very  conveniently  studied  by  the  Monte  Carlo  mi'tnod.  Indeed,  some  ol 
the  important  phenomena  can  only  be  studied  by  this  technique.  in 
addition,  the  study  of  electron  transport  in  these  types  of  structures 
requires  a  more  complete  knowledge  of  the  sub-micron  related  phenomena 
mentioned  previously. 

A  major  goal  of  this  research  program  is  to  examine  the  Monte 
Carlo  method  and  its  general  applicability  to  the  study  of  small  device 
physics  and  phenomena.  It  is  strongly  believed  that  of  the  possible 
techniques  to  be  appropriate  for  such  studies  the  Monte  Carlo  method 
will  again,  as  in  tlie  past,  prove  to  be  the  more  complete  and  accurate 
computational  method  to  study  transport  and  hot  electron  phenomena  in  the 
next  generation  of  electron  devices. 


2.0  KKSKAKCH  RESULTS 


2.1  Rcvifw  i>l  I’rogri'ss 

This  siitiiMi  I'oiuaiiis  .1  bricl  ci  i  sciiss  li>ii  ol  tliii  work  per  1  oriiifd 
during  the  [)ast  ve.ir.  .Some  ot  this  work  represeiiL.s  the  conclusion 
of  efforts  begun  more  than  .1  ye.ir  ago,  and  tliese  are  piesenled  here  lor 
completeness.  Also,  a  discussion  ot  recently  initiated  work  is  Included. 
During  the  past  year  we  have  attacked  some  problems  which  were  felt  as 
barriers  toward  the  development  of  the  Monte  Carlo  method  for  the  study 
of  small  device  effects.  Included  in  this  section  is  a  list  of 
publications  and  presentations.  Finally,  the  appendices  Include 
manuscripts  published  during  the  year. 


1)  The  work  on  central  1(000)  valley  negative  resistance  has  been 
completed  and  published  recently  in  Solid-State  Electronics  (Appendix  4.1). 
This  work  is  mentioned  here  because  it  has  basic  importance  for  velocity 
overshoot.  We  have  studied  velocity  overshoot  in  materials  which  are 
dominated  by  central  valley  transport,  and  a  paper  will  be  prepared  for 
publication.  This  particular  phenomena  could  be  a  useful  effect  for 
devices.  However,  the  materials  growth  techniques  for  the  materials 
probably  are  not  adequately  developed,  since,  for  example,  y^^y 

is  a  material  which  is  dominated  to  a  large  extent  by  central-valley 
transport.  The  phenomena  cannot  be  studied  in  the  quaternary  at  this 


time  because  of  materials  growtli  problems. 
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2)  We  have  examined  carefully  and  critically  the  definitions  ot 
estimators  for  the  diffusivity  at  high  electric  fields  in  GaAs .  The 
reasons  for  this  are  due,  first,  to  tlie  very  large  disagreement  between 
tlie  measured  diffusivity  tor  electrons  ir  GaAs  by  the  t ime-of-f light 
techniciue  (19  1  ‘‘ud  theoretical  Monte  Garlo  calculations  [20].  This 
basic  disagreement  has  not  been  satisfactorily  answered  (e.g.  see 
references  in  Appendix  A. 2).  Secondly,  it  appears  that  as  one  goes  to 
the  snui  1 1  device  dimensions  diffusion  will  become  more  and  more  Important 
as  a  possible  physical  limitation  to  device  performance.  Thus,  we  felt 
Chat  it  was  important  Co  establish  once  and  for  all  the  validity  of  the 
Monte  Carlo  method  for  Che  calculation  of  high-field  diffusion  coefficient 
in  compound  semiconductors.  This  has  been  accomplished,  since  it  has  been 
concluded  that  the  large  discrepancy  between  experiment  and  theory  is 

due  to  unaccounted  for  circuit  effects  in  the  original  experiment.  A 
paper  on  this  problem  has  been  accepted  for  publication  in  Solid  State 
Electronics  (see  Appendix  A. 2). 

3)  The  theoretical  modeling  of  alloy  scattering  has  been  continued, 

and  another  paper  on  this  subject  has  been  published  recently  (Appendix  A. 3). 
At  the  present  time  we  can  only  conclude  that  an  alloy  scattering  mechanism, 
or  some  other  form  of  scattering  similar  to  alloy  scattering  is  essentially 
dominating  the  transport  in  the  quaternary,  and  we  are  investigating  other 
mechanisms  for  their  effects  [13].  This  work  will  continue  since  we  have 
an  experimental  program  dealing  with  the  growth  and  characterization  of 
the  quaternary,  and  the  data  obtained  in  this  program  will  be  examined 
using  our  theoretical  models. 
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4)  We  have  also  been  exaniinln^;  the  eltevtive  mass  interpolation  in  the 

Ga,  In  F,  As  system,  in  an  effort  to  explain  the  contradictions  with 
1-x  X  1-y  y  ' 

experimental  effective  mass  values  obtained  by  Shubnikov-de  Hass  oscilla¬ 
tions  (see  Introduction).  We  are  presently  exploring  the  possibility  ttiat 
a  b.isic  error  has  been  matle  in  extracting  effective  mass  values  from  the 
experimental  data.  Our  initial  results  show  that  the  correction  required 
to  extr.ict  the  h.ind-edge  ellective  mass  from  the  mass  measured  at  the 
Fermi  energy  at  4.2  K  is  much  larger  than  reported,  'ihus,  the  band-edge 
corrected  etfective  masses  are  much  smaller  than  originally  reported  and 
are  much  closer  to  the  interpolated  values  than  the  experiments  have 
originally  indicated.  We  have  collaborated  with  Drs.  J.  B.  Restorff  and 
Bland  Houston  at  the  Naval  Surface  Weapons  Center  on  this  problem,  and 
we  plan  to  submit  a  publication  in  the  near  future. 

5)  The  work  begun  almost  two  years  ago  on  the  inclusion  of  the  Pauli 
Exclusion  Principle  in  the  Monte  Carlo  method  is  nearing  completion, 

and  a  paper  is  in  the  process  of  being  completed  and  submitted  for  publi¬ 
cation  (see  Appendix  4.4).  Tlie  method  is  somewhat  more  general  than 
previous  attempts,  and  has  been  used  to  calculate  low-field  drift  mobility, 
low-field  Hall  mobility,  and  high-field  drift-velocity.  We  view  this  work 
as  the  first  step  toward  including  electron-electron  interactions  as  a 
basic  mechanism  in  the  Monte  Carlo  method. 

6)  The  Monte  Carlo  program  has  been  modified  to  allow  the  addition 

of  a  sinusoidal  (in  time)  electric  field  and  the  calculation  of  mobility 


(amplitude  and  phase)  as  a  function  of  frequency.  Some  preliminary 


results  have  been  obtained  for  GaAs .  For  a  dc  (bias)  field  E 


728A  obtained  from  tlie  standard  (dc) 


These  data  are  described  reasonably  well  by  a  simple  single  pole  model 


with  a  corner  (-3dB)  frequency  of  100  GHz 


We  shall  continue  these  studies  and  publish  the  results  during  the 


coming  year 


7)  The  investigation  of  transient  (step  field)  response  is  continuing 


An  "ensemble"  Monte  Carlo  program  is  used  to  compile  carrier  histories 


(position,  energy,  etc.)  from  which  the  spatial  distribution,  k-space 


distribution,  and  various  instantaneous  ensemble  averages  (velocity 


energy,  position,  variance  of  positions,  etc.)  can  be  determined.  Results 


have  been  obtained  for  n-GaAs  (N 


and  20  kV/cm  representing  the  ohmic,  threshold,  negative  resistance,  and 


10 

A  paper  describing  these  results  is  in  preparation.  In  addition,  we  liave 
examined  tlie  influence  of  effective  mass  on  velocity  overshoot  in 
Ca^  x'''x''l  have  puljlished  a  paper  on  tliis  work  (Appendix  4.5). 

file  conclusions  ol  this  [)aper  deiiuuis t rat e  tlie  real  need  to  resolve  the 
question  t)l  the  validity  ol  the  i  liter  i)u  lat  ion  procedure  for  predicting 
effective  mass  in  the  quaternary. 


8)  We  have  initiated  a  study  ol  the  phenomena  of  real-space  electron 
transfer,  which  is  a  concept  originated  by  Hess  and  Streetman  at  the 
University  of  Illinois  [21].  The  work  performed  here  has  been  done  in 
collaboration  with  this  research  group  and  is  proving  to  be  a  very 
fruitful  interaction.  Real  space  electron  transfer  involves  the  use  of 
layered  heterostructures  of  the  type  recently  used  to  study  modulation 
doping  in  superlattice-type  materials  [22].  The  structures  employ  a  set 
of  alternating  lattice-matched  semiconductor  layers  with  dissimilar 
energy  bandgaps  and  electron  mobilities,  as  schematically  depicted  in 
Figure  1  using  GaAs  and  Ga^^  ^Al^As,  For  an  applied  electric  field  parallel 
to  the  interface,  the  electrons  are  initially  primarily  located  in  the 
high  mobility  material.  As  the  field  increases  the  electrons  can  be 
thermionically  emitted  into  the  low  mobility  material  when  their  average 
energy  exceeas  the  heterojunction  energy  barrier.  This  represents  a 
mechanism  for  negative  differential  mobility  which  is  similar  to  the 
Gunn  effect  mechanism  but  which  occurs  in  real  space  rather  than  in  k-space. 
Figure  2  shows  the  results  of  a  preliminary  calculation  for  real  space  electron 
transfer  which  illustrates  that  the  Monte  Carlo  method  predicts  the  process 
to  be  a  real  effect.  We  are  continuing  to  work  on  many  aspects  of  this 
problem.  Some  phases  of  this  work  will  be  funded  by  the  U.S.  Army  Research 
Of  f I cc , 


9)  Work  on  the  two-Jimen.s  ionul  modeling  ol  I'KT  deviees  Is  continuing. 

A  general  two-dimensional  device  model  has  been  developed,  and  during 
the  past  year  a  considerable  el  fort  has  been  expended  in  studying  the  use 
of  finite  difference  numerical  algorithms  to  solve  the  semiconductor 
transport  equations.  We  are  |)resently  not  satisfied  with  the  results 
achievable  from  this  teclinlque,  primarily  regarding  computational 
accuracy  and  convergence.  A  technique  has  now  been  developed  which 
combines  the  finite  difference  technique  witii  finite  element  methods. 

It  appears  that  the  finite  element  metliod  is  now  evolving  as  the  more 
appropriate  technique  for  the  modeling  of  complicated  device  structures 
[23].  Thus,  during  the  next  year  we  plan  to  move  in  this  direction. 

The  long-term  goal  of  this  work  is  to  provide  a  device  program  which 
will  interact  with  the  Monte  Carlo  program  to  provide  a  self-consistent 
solution  for  transport  and  terminal  properties  of  sub-micron  devices. 
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Grown  by  Liquid-Phase  Epitaxy",  Inst.  Phys.  Conf.  §er^  No^  4^, 
pp  239-247,  January,  1979. 

2.  J.  R.  Hauser,  T.  H.  Gllsson,  and  M.  A.  Littlejohn,  "Negative 
Resistance  and  Peak  Velocity  in  The  Central  (000)  Valley  of  III-V 
Semiconductors",  Solid  State  Electronics,  vol.  22,  pp.  487-493, 

July,  1979. 

3.  M.  A.  Littlejohn,  L.  A.  Arledge,  T.  H.  Glisson,  and  J.  R.  Hauser, 
"Influence  of  Central  Valley  Effective  Mass  and  Alloy  Scattering 

on  Transient  Drift  Velocity  in  Ga.  In  P,  As  ",  Electronics  Letters, 
vol.  15,  pp.  586-588,  September,  19^9.^  i  y  y 
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B.  Presentations  During  1979 

1.  M.  A.  Littlejohn,  J.  R.  Hauser,  and  T.  H.  Glisson,  "Ill-V  Compound 
Semiconductors  For  Future  Flectronic  Devices",  Invited  paper  presented 
at  the  Third  Biennial  Univer si ty/ Industry/Government  Microelectronics 
Symposium,  Texas  Tech  University,  Lubbock,  Texas,  May  21-23,  1979. 
Sponsored  by  NSF/IKEE. 

2.  T.  11.  (ilisson,  J.  R.  Hauser,  and  M.  A.  Littlejohn,  "Monte  Carlo 
Simulations  ol  Real  Space  Transfer",  Paper  presented  at  the  MBE 
Workshop,  University  ol  Illinois,  Urbana,  Illinois,  September  20-21, 
1979.  Sponsored  by  ONR/ARO. 
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C.  Paper  Accepted  For  Publication 

1.  T.  H.  Glisson,  R.  A.  Sadler,  J.  R.  Hauser,  and  M.  A.  Littlejohn, 
"Circuit  Effects  in  Time-of-Flight  Diffusivity  Measurements", 
accepted  for  publication  in  Solid  State  Electronics. 
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NKC.ATIVH  RHSISTANCK  AND  PEAK  VEI.OCITY  IN 
THE  CENTRAI  ((HH))  VAI  I.EY  OF  III-V  SKMK’ONDUCTORSt 

J  K  MAlishR,  T  M  (liissitN  and  M  A  l.irinjtiMN 
I  lri.liK.ll  I'nidnrriinK  Dcpatlnirnl,  N  ('  Stale  llniversily.  Kr'-ifih.  NC  ’7NI7,  USA 

iHi'iftvrJ  7  Junr  IS7H.  m  rrrnfil  fttrm  M)  Soirmhff  IV7KI 

Abilrarl  The  nriialivc  rrsislancr  in  III  V  malriials  such  as  (iaA^  at  large  clrrlrK  ticld\  is  generally  rerogni/eil  as 
arising  from  the  transfer  of  elei  Irons  from  the  senlral  (000)  valley  to  highei  lying  minima  in  the  conduction  hand 
Monte  C  arlo  transport  studies  show  that  the  negative  resistance  effect  is  still  present  in  III-V  materials  when  the 
valley  spacing  is  increased  to  large  values  (  'O.CeV)  and  even  present  when  the  higher  minima  are  eliminated 
entirely  from  the  calculations  This  negative  resistance  arises  from  basic  transport  properties  of  the  central  valley 
Studies  are  presented  of  the  basic  negative  resistance  effect  in  the  central  valley  of  III-V  materials  as  well  as 
studies  of  All  .In.. As  Is  -■0  7C|  and  Cia.  .In.  As  |i  O.h)  which  ate  two  specific  materials  where  the  negative 
resistance  effect  is  due  predcminantly  to  the  central  valley 


1.  ivruoowTioA 

The  negative  resistance  effect  as  first  reported  by 
Gunn|l.2|  has  been  observed  in  several  of  the  III-V 
binary  and  ternary  compounds  that  are  direct  bandgap 
materials.  The  major  factors  determining  the  negative 
resistance  characteristic  are  now  thought  to  be  fairly  well 
undershHxj.  and  a  good  summary  has  recently  been  given 
by  Ridley |3|.  A  brief  summary  of  the  conventional 
understanding  is  useful  in  relationship  to  the  present 
work. 

The  negative  resistance  effect  in  GaAs  has  been 
identified  with  a  transfer  of  electrons  at  large  electric 
fields  from  the  high  mobility  T  or  central  valley  to  higher 
lying  L  or  X  minima  having  lower  mobilities.  The  vali¬ 
dity  of  this  basic  model  has  been  verified  by  simple 
analytical  calculations|.3-5|  as  well  as  by  more  detailed 
Monte  Carlo  calculations  |fi- II).  The  energy  bands  for  a 
general  III-V  semiconductor  are  shown  in  Fig.  I.  where 
Af  is  the  energy  separation  between  the  central  valley 
and  the  next  lowest  valley,  which  is  assumed  to  be  the  /, 
valley.  In  some  materials  a  three  band  model  including 
both  the  .Y  and  /.  minima  is  required  to  give  giHid 
agreement  between  theory  and  experimental  results. 

In  the  presence  of  a  large  electric  field,  central  valley 
electrons  are  heated  by  the  field  to  large  kinetic  energies, 
and  as  the  average  kinetic  energy  approaches  the  valley 
separation,  a  large  percentage  of  the  electrons  transfer  to 
the  upper  valley.  This  transfer  is  enhanced  by  the  larger 
effective  mass  of  the  upper  valleys  as  compared  with  the 
central  valley.  Ridley  has  used  an  energy  balance  ap¬ 
proach  to  estimate  the  average  energy  of  the  electrons  as 
a  function  of  electric  field  and  equated  this  average 
energy  to  the  valley  separation  to  estimate  the  threshold 
electric  field  for  a  variety  of  transferred  electron 
materials|.1| 

Simple  considerations  of  the  transferred  electron  effect 


♦This  work  was  supported  by  a  research  grant  from  the  Office 
of  Naval  Research.  Arlington.  Virginia 


Fig.  I.  Energy  band  diagram  for  general  III-V  direct  bandgap 
semiconductor 

lead  to  the  conclusion  that  high  peak  velocities  should 
be  achieved  in  materials  with  >a)  a  large  low  field  mobil¬ 
ity  and  (b)  a  large  energy  separation  between  the  central 
valley  and  the  higher  lying  minima.  The  binary  semicon¬ 
ductor  GaAs  has  a  relatively  small  valley  separation  on 
the  order  of  0.31  eV.  The  separation  for  InP  is  somewhat 
larger  at  0.60  eV  and  InP  is  predicted  to  have  a  larger 
peak  velocity  than  GaAs  although  the  field  at  which  the 
peak  velocity  occurs  is  also  somewhat  larger. 

The  largest  valley  separation  in  the  binary  III-V 
materials  is  l.ll  eV  in  In  As.  This  large  valley  separation, 
however,  cannot  be  used  in  transferred  electron  devices 
fTED's)  because  of  the  low  breakdown  voltage  of  InAs 
due  to  the  small  bandgap  of  0..36eV.  Thus,  of  the  binary 
III-V  materials,  InP  is  predicted  to  have  the  largest  peak 
velocity  of  the  materials  which  are  useful  for  TEDs. 

The  ternary  and  quaternary  III-V  materials  provide  a 
wider  range  of  bandgaps  and  valley  separations  for 
potentially  higher  peak  veliKities  than  can  be  achieved  in 
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lh<‘  hiiiarv  itulcruiU  I  he  qii.iletnaiv  (i.i,  .In.P,  ,  Xs. 
lallKC  makhcil  li>  Ini’  Mih\li.ilc\  has  tiM'iilIv  tx-eii  pic 
JkIcO  111  have  a  peal  veliKilv  lai|:ei  Ih.iii  lhal  ol  eilhei 
(laAv  Of  I  he  lai|;e  lou  h'  lil  iiiohililv  iieetleil  toi 

hiKh  veliKiliev  leipiiiev  low  hamlfiap  iiialen.iK  with  small 
etfcelive  masses  Howevei  Iihi  low  a  hamleap.  .is  m 
IriAs.  le.iiJs  lo  low  hie.ikilown  voll.iees  thus  .1  lom 
piomise  imisl  tv  foiiiul  Ivlweeii  these  Iwo  leiimiemenis 
f  rom  hreakilowii  voll.i^e  .mil  ileviie  i oiisuh'i alions.  ihe 
h.iiHtii.ip  Ilf  hi^h  veltKilv  maleii.ils  loi  eilhei  II  IK  01 
I  K  I  s  shoulil  prohahiv  tv  .iioiiiul  I  e\  01  l.iiuei 

If  a  minimum  haiidnap  of  aioiiiul  I  e\  is  seliMeil  then 
all  Ihe  lll-V  teriiaiv  ami  ijualeiiiaiv  m.ileii.ils  i.iii  tv 
searehcii  lo  find  which  m.ilerial  has  the  laitiesi  valley 
separation  This  malenal  turns  out  lo  tv  Ali  .In,  As  for 
vyhich  Ihe  energy  hand  di.iyir.im  vs  composilion  is  shown 
in  Fig  2  .At  a  composilion  of  1  0  the  energy  gap  is 

about  091  eV  while  Ihe  F  lo  /  valley  separation  is  about 
l.l2eV.  Monie  Carlo  calculations  for  this  composition 
(to  be  discussed  in  the  next  section)  lead  to  a  low  held 
mobility  of  about  1 1.700cm‘7V  sec  and  this  coupled  with 
the  large  valley  separation  leads  lo  an  expected  l.irge 
peak  vekvily  before  Ihe  onset  of  a  negative  resistance 
effect. 

Another  ternary  material  which  has  a  large  valley 
separation  and  also  has  a  large  bandgap  is  (la,  ,ln,As 
with  X  O.A-O.h.  l  or  example  at  v  0  4  the  bandgap  is 
OSbeV  and  Ihe  valley  separation  is  0  72eV  These 
values  arc  not  quite  as  giHHl  as  those  for  Ihe  Al,  .In,. As 
system  but  Ihe  technology  for  (la.  .In,. As  is  much  more 
highly  developed  The  quaternary  (la,  .In. I’,  .As, 
previously  reported  on  1 121  also  has  a  favorable  bandg.ip 
and  valley  separation  for  l.irge  peak  velocities 

Monte  Carlo  calculations  on  both  Ihe  AT  .In,  As  and 
Ga,  .In, As  systems  have  not  shown  as  large  a  (vak 
velivily  as  initially  expected.  This  has  been  found  lo  be 
due  to  fundamental  physical  limitations  for  Ihe  peak 
velocity  which  cveur  within  Ihe  central  valley  .  For  valley 
separations  larger  than  about  0.,^eV,  Ihe  peak  velocity 
and  the  threshold  field  have  been  found  lo  be  determined 
almost  entirely  by  Ihe  properties  of  the  central  valley. 
Materials  with  large  valley  separations  will  he  referred  lo 
in  this  work  as  central  valley  dominated  materials  since 
such  materials  show  a  peak  velocity  and  a  negative 
resistance  determined  mainly  by  the  central  valley.  The 
high  field  properties  of  such  materials  arc  discussed  in 
detail  in  the  next  section 

2.  CgNTBAI.  VAI  I.HV  nrlMINATgl) 

MATyJIIAl.S 

The  calculated  velocity  field  curve  for  .AIo!iln,i7iAs  is 
shown  in  Fig.  .V  The  material  parameters  used  in  Ihe 


tThe  possibility  exists  in  the  ternary  materials  of  two  l.O 
phonon  modes,  while  in  the  present  work  a  single  I  (1  phonon 
energy  is  used  This  has  been  estimated  from  the  I  O  phonon 
energies  of  the  binary  materials  and  should  be  considered  as 
some  average  phonon  energy  Detailed  properties  of  Ihe  ternary 
materials  are  not  sufficiently  known  al  present  lo  iiislify  Ihe  use 
of  two  1.0  phonons  However.  Ihe  use  of  a  single  1  (t  phonon 
should  not  detract  from  Ihe  mam  conclusions  of  this  work 


Tig  2  Calculated  energy  gaps  vs  composition  for  All  .In,  As 

calculations  were  obtained  from  the  binary  lll-V 
parameters  as  discussed  in  previous  publications!  10- 12] 
and  are  given  in  Table  I.I  Included  in  the  calculations 
are  scattering  processes  due  to  acoustic  phonons,  polar 
optical  phonons,  ionized  impurities  (IO"'/cm'  impurity 
density),  piezoelectric,  alloy  scattering  processes,  and 
equivalent  and  nonequivalent  intervalley  processes.  The 
calculated  peak  velocity  is  about  2.7  x  10’ cm/sec  at  a 
field  of  aKiut  4(KX)  V/cm.  This  peak  value  is  about 
higher  than  the  calculated  and  measured  values  for  GaAs 
and  this  verifies  lo  some  extent  Ihe  discussion  of  Ihe 
previous  section  as  lo  Ihe  need  for  a  large  valley  separa¬ 
tion  to  achieve  large  peak  veliKilies.  However,  Ihe  peak 
velocity  is  not  as  large  as  was  expected 
The  numerical  values  shown  at  various  points  along 
Ihe  curve  of  Fig.  A  give  the  percentages  of  electrons  in 
the  /.  plus  .V  minima.  Al  Ihe  peak  of  the  vekK'ily  curve  it 
is  seen  lhal  only  0.1%  of  the  electrons  have  been  trans¬ 
ferred  lo  the  upper  valleys.  At  l(T  V/cm  which  is  far  into 
the  negative  resistance  region  only  9.2%  of  Ihe  electrons 
have  been  transferred  to  Ihe  upper  valleys.  These  per¬ 
centages  which  were  obtained  from  the  Monte  Carlo 
calculations  are  much  too  low  to  account  for  Ihe  peak 
velocity  and  the  negative  resistance  for  this  material.  For 


tCECTRIC  FIELD  1  V/cm) 

Fig  A  t'alculaled  velocity-field  curve  (or  AlcidnonAs.  The 
values  along  Ihe  curve  show  Ihe  percentages  of  electrons  in  Ihe 
upper  (111  and  100)  valleys 
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example,  at  l(FV/cm  a  transfer  of  all  electrons  hack  to 
the  central  valley  would  increase  the  vekKity  by  no  more 
than  10%  and  this  is  much  tiH>  small  an  effect  to  eli 
minate  the  negative  resistance.  This  leads  to  the  con¬ 
clusion  that  the  peak  veliHity  is  being  controlled  by  the 
central  valley  in  this  material  and  is  not  due  to  electron 
transfer  to  the  upper  valley. 

The  dominance  of  the  central  valley  in  determining  the 
peak  velocity  is  shown  in  Fig.  4  which  compares  the 
calculated  velocity-held  relationship  of  Fig.  3  with  that 
obtained  by  including  only  the  central  valley  in  the 
Monte  Carlo  calculations.  The  two  calculations  give  very 
similar  velocity-field  relationships  with  the  peak  velocity 
being  unchanged  by  the  elimination  of  scattering  to  the 
upper  valleys,  A  comparison  of  the  two  curves  in  Fig.  4 
shows  that  transfers  to  the  upper  valleys  increase  the 
magnitude  of  the  negative  resistance  beyond  the  peak 
but  the  upper  valleys  are  not  required  for  the  existance 
of  the  negative  resistance.  Upper  valley  transfers  are 
also  seen  to  cause  a  larger  drift  velocity  at  very  large 
fields  (at  1(3’  V/cm  for  example). 

A  similar  negative  resistance  effect  has  been  predicted 
to  occur  in  the  central  valley  of  InSb  at  low  temperatures 
(7T’K)(I3|.  However,  it  does  not  appear  to  have  been 
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previously  recognized  that  the  negative  resistance  effect 
of  certain  ternary  materials  can  be  dominated  by  the 
central  valley  even  at  riHtm  temperature. 

The  dominance  of  the  central  valley  is  not  limited  to 
just  .M,  .In. As  but  has  also  been  seen  in  other  ternary 
materials  such  as  Oai  .In, As  with  .i'0.4-0.6.  The 


10*  I0>  10*  10* 

ELECTHIC  field  (V/cml 


Fig.  4  Comparison  of  velocity-field  curves  for  Aloi.lnonAs 
using  central  valley  only  and  using  all  valleys. 
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vrliKilv  ficlJ  turvc  fot  ihis  Icin.irv  wilh  i  ((h  iv  slumii 
in  H(!  '  This  p.irlicul.ir  ni.ilrii.il  h.i\  .in  riiri>;v  ti.ip  of 
Oh^cV  and  .1  v.illr^  srp.ir.iln'ii  of  0  4|lrS  As  *ilh 
Alo  As.  Ihr  presence  of  (he  iip[H-i  c.illevs  h.is 

essenlullv  no  effect  on  the  pe.ik  celocilc  or  llireshoM 
field  for  netiafne  resistance  The  iippei  c.ille\s  .le.iiii 
increase  the  neiialise  resisfance  elfesi  and  iiicie.ise  itie 
hiifh  field  veliKifv 

Ihe  two  m.ilerials  discussed  so  far  hace  laiKe  s.illev 
separations,  .ind  this  is  Ihe  coiidilioii  tor  Ihe  (H'.ik  selo 
cilv  and  Ihreshold  held  beini:  deleiniined  by  Ihe  ceiilral 
valley  lor  lower  valley  sep.ir.ilioiis.  however,  such  .is 
(he  (I'leV  fvir  (i.iAs  ihe  mieslion  arises  as  lo  how 
iniporlani  is  Ihe  inlerv.illev  li.insfer  of  eleciroiis  in 
delerminini:  Ihe  pe.ik  velosilv  (  .ilciilalions  for  Ci.iAs 
have  shown  lhal  wilhoiil  the  /  .iiid  \  v.illeys  pieseiil  Ihe 
velocity  would  peak  .il  .iNuil  ’  '  s  ||i  vui/sec  .is  opposed 
lo  the  value  of  aNiiit  ’Os  10  cni/sec  including;  Ihe  iippei 
valleys  Thus  the  iransfer  of  c.irriers  lo  Ihe  upiHu  v.illevs 
reduces  Ihe  peak  velocilv  in  (ia.As  by  aboul  ’O'; 

[he  central  valley  donunaled  materials  such  as  shown 
in  h'lg  4  are  predicted  by  the  Monte  Carlo  calculations  to 
have  a  large  peak  to-valley  ratio  for  the  velocity  field 
curve  The  peak-lo-vallc)  ratio  in  l  ig  4  ranges  from 
about  4  to  7  for  peak  fields  from  20  lo  100  k\  For  (Ia.As 
with  a  smaller  valley  spacing.  Ihe  peak  to  valley  ratio  is 
at  most  aboul  2.'  by  both  experimental  measurement  1 14| 
and  theoretical  calculalionsllil  This  large  peak-lo  valley 
ratio  for  the  central  valley  dominated  semiconduclors 
may  be  especially  useful  in  certain  device  applicalions 

The  preceding  discussion  has  demonstrated  lhal  when 
the  valley  spacing  becomes  large  the  peak  vclvKity  and 
threshold  field  are  no  longer  determined  by  electron 
transfer  from  the  central  valley  lo  higher  lying  minima. 
The  peak  veliKily  is  rather  determined  by  Ihe  fun 
damental  transport  properties  of  ihe  central  valley  In 
order  lo  understand  these  central  valley  limitations  a 
study  has  been  made  of  the  high  field  properties  of  moilel 
semiconduclors  with  only  a  central  valley  These  studies 
are  discussed  in  detail  in  the  next  section 

1  riXX'VS«*T  IN  ONfRAI  yxi,i.KV  ONiy 

This  section  consulers  a  semiconduclor  in  which  the 
spacing  between  ihe  central  valley  |T((XX)t  valley!  and 


higher  lying  iiiinuiia  (either  V  or  /  valley)  is  so  large  lhal 
oiilv  Ihe  central  v.dley  must  be  considered  in  Ihe  (rails 
poll  In  all  of  (lie  III  leriiary  and  qu.ilernary  iii.ileri.ils 
Ihe  spamig  s.iii  never  be  made  so  l.irge  lhal  il  has  no 
etfeci  on  Ihe  lianspoil  process  However,  as  Ihe 
previous  section  ilhisirales,  ihe  spacing  can  be  so  large 
lhal  iiileiv.illey  liansfeis  have  a  negligible  etfeci  on  Ihe 
IH'.ik  velocity  and  Ihreshold  field  for  ccriain  materials  .A 
study  of  iraiisporl  in  only  Ihe  central  valley  can  thus  lead 
lo  an  u lute r si. Hiding  of  Ihe  physical  processes  leading  lo 
a  peak  velocity  and  negative  lesistance  etfeci  in  these 
malerials  as  well  as  provide  .111  upper  limit  lo  Ihe  peak 
veUicily  foi  olhei  materials  Since  Ihe  AI„.slnojsAs 
comes  closesi  of  .ill  Ihe  leriiary  III  V  malerials  lo  being 
domuialed  by  Ihe  central  valley.  Ihe  parameters  for  the 
c.ilculalioiis  pteseiiled  in  this  section  using  only  Ihe  cen 
Ir.il  v.illev  .lie  selected  as  those  of  this  material,  and  are 
listed  III  fable  I 

The  negative  lesisl.ince  and  peak  velocity  in  the  central 
v.illey  have  been  found  lo  be  due  lo  Ihe  fundamental 
properties  of  Ihe  nonparatxilic  bands  of  Ihe  central  valley 
as  first  discussed  by  Persky  and  Bartelinkl Ifl.  The 
negative  resistance  occurs  when  the  nonparabolic  bands 
.ire  considered  with  only  Ihe  dominant  scattering  process 
of  polar  optical  scattering.  This  is  seen  in  F'ig.  6  which 
compares  Ihe  velocity-field  curve  for  only  polar  optical 
scattering  with  the  corresponding  curve  including  all 
scattering  prtKesses.  The  general  shape  of  the  curves  is 
seen  lo  be  determined  by  polar  optical  scattering  with  Ihe 
other  scattering  processes  simply  reducing  the  velocity 
by  almost  a  constant  factor  at  all  field  values.  This 
dominance  of  polar  optical  scattering  is  to  be  expected 
since  it  is  the  major  energy  loss  process. 

The  presence  of  a  negative  resistance  in  Monte  Carlo 
calculations  when  only  polar  optical  scattering  is  con¬ 
sidered  is  at  first  very  surprising.  This  differs  from  pre¬ 
vious  studies  where  it  has  been  shown  that  only  polar 
optical  scattering  leads  to  an  energy  runaway  condition 
at  large  electric  fields|l6-l8|.  .As  shown  below  the 
difference  between  these  earlier  studies  and  Ihe  present 
study  arises  from  Ihe  description  of  Ihe  energy  bands. 
The  studies  referenced  above  considered  parabolic 
energy  bands  while  Ihe  present  work  has  considered  a 
non  parabolic  band  described  by 


Fig  .V  VeliKitv  field  curves  for  (iaosInojAs 


r — if  =  r(e)  =  e(l -foe).  (I) 

2m* 

where  m*  is  the  effective  mass  at  e  =  0  and  a  is  the  first 
order  non-parabolicity  facIor|l9|.  With  parabolic  bands 
one  gets  energy  runaway  and  with  non-parabolic  bands 
one  does  not. 

The  basic  physics  of  high  field  transport  can  be  des¬ 
cribed  analytically  by  energy  and  momentum  balance 
equations  of  (.1! 
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where  h  is  Ihe  eleclrk  held,  m*  the  eneruy  dependent 
etfeelive  mass.  »,  ihe  Ihernial  cneriiy.  Ihe  nuinienUini 
relaxalion  lime  and  r.  Ihe  enerKV  rehualinn  lime  Ihese 
ean  he  eomhined  .ind  soliilinns  uhlained  fur  /■  and  r  as 
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If  a  model  is  available  for  t„  and  r,  as  funelions  of 
energy,  these  edualions  express  /•'  and  r  as  funelions  of 
average  electron  energy.  Hy  solving  these  for  various 
values  of  energy,  e  as  a  function  of  f'  can  be  obtained 
Fur  polar  optical  scattering  only  and  at  large  energies 
|(  *■  Aiuii)  the  relaxation  rates  can  he  approximated  as|^|. 
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f  iR  VeliKify  field  curve  for  cenlrul  valley  wilh  just  fH>lar 
optical  scattering. 


t  -fj  _  (/^2m*)''^(h<i>o)’  / J _ [^ \  I  d> 

T,  SttA"  \«.  fu/y(7?^df  (6) 

xln[4y(e)(hn*,‘^’)  ’  j,  (7) 

where  fhuo  is  the  optical  phonon  energy  and  n(iao)  is  the 
Bose-Einstein  factor.  More  detailed  and  complete  treat¬ 
ments  of  the  relaxation  rates  can  be  found  in  Refs.  1 20. 2 1 1. 
However,  Ihe  present  simple  expressions  appear  to  contain 
all  the  information  needed  to  understand  the  central  valley 
negative  resistance  The  energy  dependent  effectiye  mass 
can  also  be  approximated  as 


In  a  parabolic  band  (y  =  r,  (dy/dt)  I)  both  Ihe  rales  of 
momentum  loss  and  energy  loss  decrease  with  energy. 
This  makes  it  impossible  to  achieve  steady  state  above  a 
certain  energy.  For  a  non-parabolic  band  as  described  by 
eqn  (I).  r„  '  approaches  a  constant  at  large  energies 
while  increases  slowly  with  energy  as 

ln(2(/A<uo).  In  this  type  of  energy  band  a  steady  state 
solution  can  be  found  at  any  energy 

The  analytical  model  as  expressed  by  eqns  (4)-(8)  has 
been  solved  for  Ihe  parameters  given  in  Table  I  and  for 
the  energy  band  of  eqn  (1).  The  results  for  velocity  as  a 
function  of  field  are  plotted  as  psrints  in  Fig.  6.  Excellent 
agreement  is  seen  between  the  results  of  the  analytical 
model  and  the  Monte  Carlo  calculations.  This  supports 
the  accuracy  of  the  Monte  Carlo  results  and  demon¬ 
strates  that  negative  resistance  is  a  feature  of  the  central 
valley  when  non-parabolic  energy  bands  are  considered. 

Two  basic  questions  can  still  be  raised  with  respect  to 
the  results  reported  here.  First,  the  energy  band  is  in 
reality  more  complicated  than  the  form  expressed  in  eqn 
(1).  This  equation  approaches  a  linear  E  vs  k  relationship 
at  large  energies  while  real  bands  have  a  region  of 
negative  effective  mass  at  large  energies  In  such  an 
energy  band,  the  electron  energy  for  a  given  field  will  he 


even  less  than  that  given  by  a  band  described  by  eqn  (I). 
Thus  if  energy  runaway  does  not  occur  with  a  band 
described  by  eqn  (I),  it  will  not  occur  in  a  band  where 
energy  increases  more  slowly  with  k.  Also  the  negative 
effective  mass  region  of  a  band  can  by  itself  lead  to  a 
negative  resistance  at  high  fields.  Thus  a  more  accurate 
description  of  the  nonparabolic  band  should  enhance  the 
negative  resistance  calculated  in  the  present  work. 

A  second  question  arises  as  to  the  neglect  of  p -state 
mixing  in  the  model  for  Ihe  relaxation  rates  of  eqns  (6) 
and  (7)|7|.  In  order  to  compare  the  Monte  Carlo  and 
analytical  results,  Ihe  Monte  Carlo  results  reported  in 
Fig.  ft  did  not  include  p -state  mixing  in  the  calculations. 
A  comparison  of  the  Monte  Carlo  results  both  including 
and  neglecting  p-state  mixing  is  shown  in  Fig.  7.  As  can 
be  seen  there  is  little  difference  in  the  overall  results.  The 
velocity  including  p -stale  mixing  is  just  slightly  larger 
than  that  neglecting  p -stale  mixing.  Thus  Ihe  existence  of 
Ihe  negative  resistance  region  is  not  sensitive  to  Ihe 
exact  model  used  for  polar  optical  scattering. 

■A  simple  interpretation  of  Ihe  positive  and  negative 
resistance  regions  can  be  obtained  from  the  momentum 
and  energy  balance  equations.  The  energy  balance  equa- 


F'ig  7  Comparison  of  calculated  velocity  field  curves  including 
iimi  not  mi'tmlinft  p-stafe  mixing  in  the  scattering  rates. 
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whfff  cgn  (I’)  i\  onlv  valid  al  largo  energies  Al  Um 
fields  and  U'w  parliele  energies  where  mj  m*  and 
arc  eonslant.  eqn  III)  prediets  a  linearly  inereasing 
veliKily  with  field  Al  large  fields  if  t„  were  eonslanl  eqn 
(12)  prediels  a  decreasing  veliKily  wilh  mcieasing  field 
The  Memle  I'atUi  results  aiul  eqn  (’)  aetually  shuw  that 
To  '  increases  with  panicle  energy  as  ln(»).  and  »  in  lurn 
increases  as  {■  '  al  large  fields  These  combine  lo  make 
ro  '  increase  as  In  If  )  This  is  a  we.ik  dependence  and 
does  nol  compensate  (or  the  (  '  slependence  of  eqn 
1 12)  Tqualion  1 1 1)  is  also  satisfied  al  large  fields  because 
m‘  increases  with  energy  and  field 
Kqualions  (II)  and  (12)  can  be  used  lo  proyide  an 
upper  limit  on  the  yeliKitv  which  can  be  achieved  in  a 
semiconductor  By  equaling  the  limiting  espressions  of 
eqns  ill)  and  (12)  and  taking  r„  r,,*  and  solying  for 
qF.rJm*.  one  obtains 


ipSi™..  =  [^"lanh(W2kr)j''‘  (1.1) 

Values  of  r™,  for  seseral  binary  and  compound  materi¬ 
als  are  shown  in  Fig  8  Also  shown  for  seseral  materials 
are  Ihe  peak  velocity  salues  that  result  from  Monte  C'iu-|o 
calculations  including  all  bands  and  all  scattering  nro- 
cesses  For  the  binary  materials  GaAs,  InP  and  InAs  the 
Monte  CiU’lo  results  give  peak  velocities  that  are  about 
82-W>  of  the  values  predicted  by  eqn  (1.1).  For  the  two 
lern-ary  and  one  qu.ilernary  materials  shown.  Ihe  Monte 
l';ulo  results  differ  from  eqn  (11)  by  large  amounts  The 


tr.M  and  ro  are  only  equal  for  special  scattering  priKCsses  such  as 
isotropic  scattering  TTiev  are  nol  exactly  equal  for  polar  optical 
scattering  However,  Ihe  limiting  expression  obtained  in  Ihi' 
manner  has  been  found  lo  provide  a  good  upper  limit  lo  the  peak 
velocilv  for  a  variety  of  semiconductors 


lower  Monte  (,'arlo  resull  for  ihe  ternary  and  quaternary 
nialerials  is  due  lo  the  presence  of  alloy  scattering  in 
these  materials  which  is  nol  present  in  (he  binary 
materials  If  alloy  scattering  had  nol  been  included  all  of 
Ihe  Monte  (  arlo  results  would  have  been  around  80 
of  Ihe  liniiling  value  given  by  eqn  111) 

Monie  Carlo  results  are  nol  shown  for  In.Sb  and  Cia.Sb 
The  band  g.ip  ol  In.Sb  is  imi  low  for  this  material  lo  be 
Used  lie. II  the  ix'.ik  velocity  Foi  (iaSb  the  piesence  of 
the  I  iiiiiiiiii.i  .It  only  OOUeV  above  Ihe  I  minima  causes 
Ihe  iieg.ilive  lesisl.iiice  effect  nol  lo  occur  in  this 
iiialeiial  Ihe  d.il.i  in  Fig  8  indicales  lhal  (here  is  a 
general  trend  toward  l.iiger  (K'ak  yelix'ilies  as  Ihe 
bandgap  of  Ihe  seniicimduclor  decreases 

ITie  fundamenlal  loles  of  polar  optical  scallenng  and 
band  nonparabolic ily  in  Ihe  cenital  valley  combine  lo 
determine  the  peak  velocity  which  can  he  achieved. 
Figure  y  shows  Ihe  velocity  field  curves  for  Ihe  central 
valley  only  of  Al,.  ;sln,i  7.As  with  different  u  values, 
where  Ihe  theoretical  value  of 


gives  ll)4eV  '  for  AL,  .,ln,i  7,,As  This  corresponds  lo 
one  of  Ihe  curves  m  Fig  V  .Also  shown  are  veliKily  field 
curves  for  a  values  an  order  of  magnitude  smaller  and  an 
order  of  magnitude  larger  than  Ihe  theoretical  value.  For 
values  of  II  less  than  0  104  the  calculated  values  are 
close  lo  Ihe  (I  104  curve  of  Fig.  .Al  Ihe  theoretical  value 
of  It  Ihe  nonparabvilicily  is  seen  lo  reduce  Ihe  k'w  field 
mobility  by  abvnil  17*^.  lo  reduce  Ihe  peak  velocity 
slightly,  and  lo  increase  Ihe  threshold  field  from  around 
IkV/cni  10  around  4kV7cm  .A  very  large  nonparaboli- 
cily  III  =  10.4  eV  '),  which  fortunately  does  not  cKcur  in 
Ihe  Ill-V  materials,  is  seen  to  drastically  reduce  the  low 
field  mobility  and  ihe  peak  veliKily  The  nonp;u-abolicity 
will  be  most  important  in  materials  with  narrow 
bandgaps  and  large  valley  spacings.  since  a  is  large  for 
such  materials  and  Ihe  carriers  become  heated  to  large 
energies  before  inlervalley  transfer  occurs. 
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4  S(  MMAK)  A>0  Um  i.lSHifiS 

The  presently  accepted  theory  of  negative  resistance 
effects  in  Ill-V  semiconductors  attributes  the  effect  to 
the  transfer  of  electrons  from  a  high  mobility  central  (I  ) 
valley  to  lower  mobility,  higher  lying  energy  valleys  (I. 
or  X).  The  present  work  has  shown  that  in  materials  with 
a  large  valley  spacing  the  negative  resistance  becomes 
dominated  by  the  central  valley  Monte  Carlo  cal¬ 
culations  have  shown  that  the  presence  of  the  upper 
valleys  is  not  required  for  a  negative  resistance  effect  but 
that  polar  optical  scattering  acting  in  a  nonparabolic 
central  valley  alone  gives  rise  to  a  peak  veUxity  and 
negative  resistance  effect  Kor  materials  with  vallev 
spacings  of  O.^eV  or  larger  the  peak  velixily  and 
threshold  field  has  been  found  to  be  determined  almost 
entirely  by  the  central  valley  negative  resistance  effect 
Transfers  to  the  upper  valleys  influence  mainly  the 
magnitude  of  the  negative  resistance  and  the  carrier 
velocity  near  the  valley  minimum  of  the  veliKity  field 
curve. 


■IV 1 

M.iterials  which  h.ivc  huge  vallev  vtsicings  and  .ire 
dominated  bv  the  centi.d  v.dley  h.ive  large  peak  to  valley 
velixitv  ratios  as  predicted  bv  Monte  (  ailo  calculations 
Ibe  tcin.uv  maten.ils  .All  .In. Asti  (17'),  (/a,  .In. As 
(I  0 ')  and  the  quaternaiv  (iai  .In. I’,  .As,(i  0  7'). 
V  0  4)  are  typical  of  central  vallev  dominated  materials. 
Ihe  largest  peak  vehxities  which  have  so  far  been 
calculated  by  the  Monte  (  arlo  technique  have  been 
found  III  these  iialeii.ils 
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Ciri-iilt  Kfft'i'ls  In  Timf-ot -I'l  iglit  Dlffuslvlty  Mcasur  c-inents 

T.  H.  Glisson,  R.  A.  SadliT,  J.  K.  Mauser,  and  M.  A.  Littlejohn 

Department  ot  Electrical  Engineering 
North  Carolina  State  University,  Raleigh,  NC  27650  | 

Abstract 

A  monte  carlo  procedure  is  used  to  examine  certain  effects  that 
interfere  with  the  determination  of  the  dlffusivity  D  from  current 
rise  and  fall  times  in  a  time-of-f light  experiment.  It  is  found  that 
electron  bombardment  time,  sample  capacitance,  contact  resistance, 
and  oscilloscope  rise  time  can  cause  considerable  error  in  the  measure¬ 
ment  of  D,  especially  near  the  threshold  electric  field.  It  is  shown 
that  a  reasonable  accounting  for  these  effects  can  explain  much  of  the 
difference  between  experimentally-and  theoretically-determined  D-E 
characteristics  for  GaAs. 


This  work  was  supported  by  the  Office  of  Naval  Research. 


I  lU  rcxluct  Ion 


I'xti'nsivo  I'HMsiir  i-mont  s  ol  t  lie  transport  properties  of  electrons  in 
CaAs  were  reputed  by  Rucb  and  Kino  in  1968  [1).  It  has  proved  difficult 
to  bring  monte -car lo  calculations  of  the  ulffuslvlty  versus  field  (D-E) 
characteristic  into  agreement  with  the  experimental  D-K  characteristic 
in  [1]  without  sacrificing  some  agreement  elsewhere,  e.g.,  with  velocity- 
field  (v-E)  or  mobility-temperature  (p-T)  measurements.  This  state  of 
affairs  is  Illustrated  in  Kig.  1,  where  three  calculated  (monte-carlo) 

D-E  characteristics  are  superimposed  on  the  experimental  (RK)  D-E 
characteristic  from  [1].  Curve  A,  calculated  by  Abe,  et  al  [2],  is  in 
good  agreement  with  the  experimental  data  for  E  ^  2kV/cm,  but  leads 
(through  the  Einstein  relation)  to  a  low-field  mobility  that  is  about 
twice  the  accepted  value.  Curve  L  was  calculated  by  the  present  authors 
using  the  model  proposed  by  Littlejohn,  et  al  [3].  That  model  yields  calculated 
v-E  and  p-T  characteristics  that  are  in  good  agreement  with  experimental 
data,  yet  yields  a  D-E  characteristic  (curve  L  in  Fig.  1)  that  is  quite 
different  from  that  measured  by  Ruch  and  Kino.  Curve  PR  in  Fig.  1  was 
calculated  by  Pozela  and  Reklaitis  [4]  using  a  model  that  is  only  slightly 
different  from  the  one  described  in  [3].  The  present  authors  have  been 
unable  to  reproduce  this  result,  and  cannot  attribute  the  great  differences 
between  curves  L  and  PR  to  the  slight  difference  between  the  models 
described  in  [3]  and  [4]. 

The  experimental  (RK)  characteristic  itself  1b  questionable  in  two 
regards.  First,  it  rises  sharply  in  the  ohmic  region  (0  <  E  <2kV), 
whereas  D  should  be  nearly  independent  of  E  there  because  of  the  Einstein 
relation.  Second,  the  RK  characteristic  appears  (when  extrapolated) 


2 


to  approach  ~100  -  150  cni*’/s  at  larj;o  fields,  wluTcas  experimental  data 

2 

Indicate  that  1)  does  not  exceed  20  cm  /s  at  K  =  50  kV/cra  [5). 

In  Ruch  and  Kino's  paper  (1)  the  experimental  D-E  characteristic 
was  compared  to  one  attributed  to  Butcher  and  Fawcett  [6]*  iind  it  was 
suggested  that  the  considerable  differences  were  due  in  part  to  "inter- 
valley  diffusion".  However,  it  is  hard  to  see  how  this  effect  could  have 
been  excluded  from  the  monte-carlo  calculation  described  in  [6],  since  a 
two-valley  model  was  used. 

In  this  paper  we  suggest  that  some  of  the  differences  between  the 
experimental  and  calculated  D-E  characteristics  might  be  due  to  seemingly 
negligible  circuit  effects  that  interfere  with  the  measurement  scheme 
described  in  [1).  These  circuit  effects  include  the  sample  capacitance, 
contact  resistance,  imperfect  matching,  oscilloscope  lise  time,  and  the 
like. 


Theory  of  Diffuslvity  Measurement 

In  the  t ime-of-f 1 ight  experiment,  carriers  are  Injected  at  the  cathode 
and  drift  under  tiie  applied  field  to  the  anode.  Diffusion  causes  the 
spatial  distribution  of  carriers  to  spread  during  the  drift.  Ruch  and  Kino 
assumed  (see  Appendix  C  of  [1])  the  spatial  distribution  of  carriers  to  be 
gausslan  with  mean  x  =  vt  and  variance  o^  =  2Dt,  where  v  and  D  are  the 
drift  velocity  and  diffuslvity,  respectively.  Thus  the  carriers  reach 
the  anode  at  times  near  1/v,  where  I  is  the  sample  thickness.  At  that 
time  o  =v2Dll/v  .  and  the  current  falls  from  95X  to  5X  of  its  peak  value 

X  ’ 

in  a  time  given  by  =  3.28  Thus  the  diffuslvity  can  be  obtained 

♦Actually,  no  D-E  data  are  given  in  [6].  The  D-E  data  attributed  in  (1]  to  [6] 
appear  to  have  been  obtained  from  the  v-E  data  in  [6]  via  the  Einstein  relation. 


from  a  measurement  of  the  95%-5%  current  fall  time  through 


where  v  Is  obtained  from  the 


e  thickness 


The  more  elaborate  development  given  in  Appendix  C  of  [1]  leads  to 


where  t  is  the  57,-957,  current  pulse  rise  time  and  v  is  obtained  from  the 


However,  it  is  not  clear  that  diffusion  effects  have  anything  at  all  to  do 


with  the  shape  of  the  leading  edge  of  the  current  pulse;  the  inclusion  of 


T  in  Eq.  (2)  is  perhaps  best  regarded  as  an  attempt  to  compensate  for 


circuit  effects  and  the  finite  time  for  which  the  sample  is  bombarded 


(to  inject  carriers).  Ruch  and  Kino  state  correctly  that  Eq ,  (2)  compensates 


exactly  for  circuit  effects  if  the  rates  of  current  rise  and  fall  and  the 


circuit  impulse  response  are  gaussian,  but  the  assumption  of  a  gaussian  cir 


cuit  impulse  response  is  hard  to  justify.  An  exponential  circuit  impulse 


response  seems  more  likely,  in  which  case  Eq.  (2)  can  give  wrong  values  for 


Circuit  Model  and  Results 


The  experimental  apparatus  and  set-up  used  by  Ruch  and  Kino  are 


described  in  [1],  where  rise  (and  fall)  times  of  "as  much  as  0.2  ns 


were  attributed  to  sample  and  support  capacitance,  contact  resistance 


and  the  like.  The  oscilloscope  rise  time  was  0.06  ns,  and  the  sample 


was  bombarded  by  electrons  (to  inject  carriers)  for  times  on  the  order 


of  0.1  ns.  The  electron-hole  pairs  were  created  within  3  pm  of  the 


cathode . 


A 

Those  effects  and  I'tliers  (such  as  sjiace  charge)  contribute  to 
and  tp  in  Eq .  (2),  perhaps  as  much  as  or  more  than  diffusion,  and  may 
distort  the  value  of  D  thus  obtained. 

A  first-order  accounting  for  circuit  effects  can  be  ma4e  by  super¬ 
imposing  the  responses  of  a  single-pole  (R-C)  circuit  to  a  number  of 
rectangular  current  pulses,  each  of  which  represents  the  contribution  of 
a  single  electron  to  the  total  current  through  (or  voltage  across)  the 
sample.  This  is  illustrated  in  Fig.  2,  where  an  individual  current  source 
is  described  by 


k(il-x  ) 

^n^"^  '  [“(t-t^)  -  u(t-t^  -  At^)].  (3) 

n 

(In  Eq.  (3),  k  is  a  dimensioned  constant  that  accounts  for  charge  density, 

sample  cross-sectional  area,  and  electronic  charge.)  In  the  calculations 

to  be  described  here,  is  randomly  selected  from  the  Interval  {0,3pm] 

(to  account  for  penetration  of  the  bombarding  electrons),  t  is  randomly 

n 

selected  from  the  interval  [0,0.1  ns]  (to  account  for  the  time  for  which 
the  sample  is  bombarded),  and  At^  is  selected  from  the  density 


£-x  -(Jl-x  -vAt  )‘’/(4DAt  ) 

r  t  r\  ^  nn  n 

f(ll-x  ,At  )  =  >  ■■  e 


(4) 


4'nDAt‘ 


which  is  the  density  of  the  time  to  absorption  (4t^)  for  Brownian  particles 
released  (at  t=0)  a  distance  to  the  left  of  an  absorbing  boundary 

(the  anode)  [7].  In  the  Ruch-Kino  experiment  and  in  the  calculations  here, 
the  sample  thickness  I  =  0.03  cm. 

Our  procedure  was  as  follows:  (1)  The  field  dependence  of  v  and  D 
in  Eq .  (4)  was  obtained  by  a  monte-carlo  calculation  from  the  model 
described  in  [3].  This  is  equivalent  to  assuming  that  curve  L  in  Fig.  1 
is  the  "correct"  one  that  would  be  observed  in  the  absence  of  Interfering 
effects;  (2)  For  each  field  (IkV/cm  <  E  <  12  kV/cm),  one  thousand  values 


5 


each  for  x  ,  t  ,  .and  At  wore  selected  as  described  above;  (3):  The 
n  n  n 

circuit  response  (simulated  pul se)  was  constructed  from 

1000  k(Jl-x  ) 

r(t)  '  )-g(t-t  -At  )]  (5) 


where 


-t/T 

g(t)  =  (l-e  ‘']u(t)  ,  (6) 

is  the  step  response  of  the  circuit  for  time  constant  (A):  The  5%-95% 

rise  time  x  ,  95%-5%  fall  time  i  ,  and  50%-50%  pulse  duration  t_  (transit 
K  r  i 

time)  of  the  simulated  pulse  r(t)  were  determined;  and  (5):  The  drift 
velocity  V  was  calculated  as  v  “  dlffusivity  was  then  calculated 

using  Eq ,  (2) . 

The  calculation  just  described  was  performed  for  several  values  of  the 

circuit  time  constant  x  .  The  best  agreement  between  the  simulated  and 

c 

experimental  D-E  characteristic  was  obtained  for  x^  =  0,25  ns.  The  results 

are  shown  in  Fig.  3.  The  dashed  curve  is  the  experimental  data  from  (1], 

and  the  solid  circles  are  the  dlf fusivities  calculated  from  Eq.  (2)  as 

described  above.  Curve  L  is  repeated  from  Fig.  1,  and  shows  the  specified 

dif fusivities  used  in  Eq.  (A)  for  the  random  selection  of  the  individual 

carrier  transit  times.  Curve  L  and  the  v-E  data  used  in  Eq.  (A)  were 

calculated  by  the  authors  using  the  model  described  in  [3]  with  impurity 

lA  -3 

concentration  N  =  10  cm  .  This  impurity  concentration  is  reasonably 
consistent  with  the  sample  resistivities  given  in  [1], 

The  observed  D-E  characteristic  in  Fig.  3  is  in  good  qualitative 
agreement  with  the  experimental  data,  but  is  displaced  toward  higher 
fields.  This  displacement  is  due  to  the  fact  that  the  model  used  leads 


-  -V 


6 


to  a  v-E  cliaracter istlc  chat  Is  displaced  to  the  right  relative  to  that 
measured  and  used  to  calculate  D  In  [1].  The  v-E  characteristic  measured 
in  [1]  is  In  much  better  agreement  with  that  in  [6]  than  that  in  [3],  so 

I 

the  calculations  described  above  were  repe-'^ted  using  the  theo^'etlcal  v-E 

and  D-E  data  attributed  in  [1]  to  Butcher  and  Fawcett  [6].  The  best  fit 

circuit  time  constant  was  x  =  0.22  ns,  and  the  results  are  shown  in 

c 

Fig.  A.  The  simulated  D-E  characteristic  is  closer  to  the  experimental 
characteristic  than  in  Fig.  3,  mainly  because  the  peaks  more  nearly 
coincide. 


Conclusions 

The  circuit  time  constants  x  =  0.25  ns  and  0.22  ns  used  to  obtain 

c 

the  results  just  described  are  reasonably  consistent  with  the  estimates 

given  in  [1]  (0.2  ns  for  sample  &  supports,  0.06  ns  for  the  oscilloscope), 

and  lead  to  observed  (measured)  diff usivities  that  are  greatly  inflated 

2 

near  their  peak  value.  The  term  x  in  Eq .  (2)  does  not  compensate  for  the 

K 

circuit  effects,  although  it  was  found  to  almost  compensate  for  the  effects 
of  the  0.1  ns  bombardment  time.  It  was  first  thought  that  the  circuit 
effects  would  also  disturb  the  transit  time  measurement,  and  thus  the 
measurement  of  v,  because  the  leading  and  trailing  pulse  edges  have  different 
slopes  and  would  suffer  different  delays;  this  was  not  found  to  be  a  problem 
for  the  time  constants  consid.;red  here.  The  calculated  drift  velocity  was 
always  within  a  few  percent  of  that  specified  in  Eq.  (4)  for  the  selection 


process. 


7 


The  impact  of  the  circuit  effects  on  the  measurement  of  D  is  further 


illustrated  by  Fig.  5.  Here  the  specified  v  and  D  were  calculated  at 


14  -3 

E  “  4kV/cm  using  the  model  in  [3]  with  N  “  10  cm  .  Thus  the  "true' 


diffusivlty  is  300  cm  /s,  as  was  observed  for  a  circuit  tim^  constant 


of  zero.  The  s imul  a  ted di f f us  ivi ty  remains  close  to  the  actual  (specified) 


diffusivlty  for  t  <  0.05  ns,  and  rises  sharply  for  >  0.15  ns.  At 


=  0.25  ns,  the  s imulated  d 1 f f us Ivl ty  is  about  900  cm  /s,  and  corresponds 


to  the  peak  value  of  the  observed  D  in  Fig.  3. 


Equation  (2)  can  give  good  results  if  the  external  circuit  effects 


are  small  compared  to  the  effects  of  diffusion.  In  terms  of  the  model 


used  here,  this  requires  that  the  circuit  time  constant  be  small 


compared  to  the  fall  time  tj,  due  to  diffusion  alone  at  all  fields  of 


interest.  The  fall  time  Tj,  is  on  the  order  of  2o^(t.j.)/v  =  2^2Dl/v  , 


and  is  smallest  near  threshold.  Based  on  this  simple  argument,  it 


appears  Eq.  (2)  is  useable  provided 


<<Y8Di/v  . 


For  G^LAs  and  E  =  4kV/cm  we  find  vSDi/v  =  0.1  ns.  As  expected,  this  value 


is  also  roughly  that  at  which  the  curve  in  Fig.  5  begins  to  turn  upward. 


If  the  Inequality  in  Eq .  (7)  is  not  satisfied,  circuit  effects  must 


somehow  be  accounted  for.  If  the  circuit  effects  can  be  modelled,  a 


procedure  like  the  one  described  in  the  previous  section  can  be  used  to 


construct  curves  like  the  one  in  Fig.  5  for  various  "actual"  values  for  D. 


These  curves  could  then  be  used  to  correct  the  results  obtained  from 


Eq.  (2), 


I 


I 

I 

r. 
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Fig.  1 

Fig.  2 
Fig.  3 

Fig.  A 
Fig.  5 


Figure  Captions 


.  Dlffuslvity  versus  applied  electric  field  for  electrons  in 
CaAs .  RK:  experimental  data  from  [l];  A:  Monte-carlo 
calculation  from  [2];  PR:  Monte-carlo  calculation  from  [4]; 

L:  Monte-carlo  calculation  by  the  present  authors  luslng  the 

model  described  in  [3]. 

.  Circuit  model  used  in  the  calculation  of  the  simulated  dlffuslvity. 
(see  text). 

.  Simulated dlf fusivlty  versus  field  when  the  "actual"  (specified 

in  Eq.  4)  D-E^and_v-E  values  are  obtained  from  the  model  in  [3] 

(with  N  =  10  cm  ).  t  =  0.25  ns. 

c 

.  Siculated dlffuslvity  versus  field  when  the  "actual"  D-E  and 

v-E  are  obtained  from  [1],  attributed  there  to  [6].  =  0.22  ns. 

.  Sin^ulated  peak  dlffuslvity  versus  circuit  time  constant.  The 
specified  dlffuslvity  (300  cm  /s)  and  drift  velocity  in  Eq.  (4) 
were  calculated  from  the  model  in  [3]. 
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Carrier  compensation  and  alloy  scattering  in 
Gai_,In,Pi_v ASv  grown  by  liquid-phase  epitaxy 
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Abxlmcf.  (  .iMKj  (iii»bihtics  aclncvetl  cxihihikiiI.iIIv  in  (ia,  ,In,P|  yAsy  alloys  arc 
coiisulcrahly  less  than  those  pietheleii  In  Monte  Uailo  tt.insptnt  eaieulations  In  this  paper,  a 
eoinpansDn  between  the  leinperaluie  ilepeiulenee  of  eleetion  I  lall  niobihty  inea\urcd  experi- 
inenlallv  on  t  I*i  -grown  samples  am!  ealeulaleil  from  a  Mt>nte  ('arlo  method  is  presented.  The 
experimental  data  is  obtamcil  lioin  samples  which  are  n-lype  and  aie  grown  on  semi -insulating 
InP  substrates  The  theoretical  ealeul.itions  are  based  on  a  ilitfusion  eoeftieient  estimator  for 
calculation  of  Hall  mobility  An  analysis  iif  the  data  mdieaies  that  both  alhry  scattering  and 
earner  compensation  are  necessary  ti>  explain  the  observed  temperature  dependence  of  the 
Hall  mobility 


1.  Introduction 

The  111-V  quaternary  alloy  Gai  ,In,Pi  ..As,,  has  received  considerable  attention 
recently  due  to  many  desirable  electronic  and  material  properties,  and  this 
compound  has  been  used  in  a  wide  variety  of  device-related  studies  (Clawson  1978). 
In  spite  of  the  wide  variety  of  developed  applications,  the  desirable  transport 
properties  proposed  for  this  material  based  on  calculations  using  the  Monte  Carlo 
method  (Littlejohn  et  al  1977a)  have  not  been  fully  realized  (Morkoc  et  al  1978, 
Houston  etal  1978).  In  particular,  the  reported  Hall  mobilities  have  been  consider¬ 
ably  lower  than  predicted,  and  room-temperature  electron  concentrations  below 
about  1  X  !()''’ cm  '  have  been  dilficult  to  achieve.  Much  of  the  work  has  b  cn  done 
using  liquid-phase  epitaxy  (i  in  ).  However,  the  electron  concentrations  are  much 
larger  for  the  quaternary  than  those  usually  achieved  for  i  I’l  growth  of  the  binary 
compounds  GaAs  and  InP.  even  though  the  starting  constituent  source  materials 
used  to  grow  the  quaternary  are  of  the  same  purity  as  those  used  for  growth  of  the 
binaries. 

In  this  paper  a  comparison  between  experimental  Hall  mobilities  of  i.PH-grown 
quaternary  layers  on  InP  substrates  and  those  calculated  from  the  Monte  Carlo 
method  is  presented.  Based  on  this  comparison,  it  is  concluded  that  both  carrier 
compensation  as  well  as  alloy  scattering  are  contributing  to  the  low-field  Hall 
mobility  and  to  its  temperature  dependence. 

2.  Experimental  aspects 

The  samples  studied  in  this  paper  were  grown  by  liquid-phase  epitaxy  at  660  °C  in  a 
standard  horizontal  slider  system  (Phatak  and  Harrison  1978  private  com- 
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iiuiiiuMlioiil  Hk'  ;is  Blown  l;isi‘is  wi-ii’  iiiulo|)cil.  anil  Ihc  sul>stiatc  malofial  was 
scnii-insulaimB  IV  ilopi'il  Ini’  in  a  ( I  I  I IH  oi  ii'iitalion.  I'vpii'al  lavoi  I  hick  nesses  were 
aroiiiul  10  pm  The  Hall  ineasnienienls  weie  inaile  in  a  lanahle-lenipeialure 
civostal  helween  77  anil  M)0  K  A  slanilaul  V'an  ilei  I’auw  sample  Beomelry  was 
iiseil.  aiul  eonlaels  to  the  samples  weie  maile  with  an  In  Sn  alloy  Ulousum  c(  iil 
I07S), 

3.  Monle  Carlo  ealculalioiis 

The  ealeulalions  piesenieil  m  iliis  papei  weie  ohiaineil  hv  the  Monle  C'ailo  melhoil. 
usinp  a  pio)>ramme  deseriheil  previously  (1  iillejohn  <•!  til  1077a).  Of  partieular 
importance  is  the  methoil  used  to  oliiam  the  i|ualeinarv  material  parameters  for 
CfUi  ,  In,  I’l  ,  As,,  This  is  ilone  hy  an  inleipolation  procedure  which  is  based  on  an 
assumed  knowledBc  of  the  constituent  Innarv  material  paiameters  of.  in  this  case. 
GaAs.  InAs.  Inl’  and  (!al*.  While  much  needs  to  he  done  to  validate  further  this 
interpolation  procedure,  it  has  proved  useful  m  calculatiiiB  eneiBV  band  jjap  and 
lattice  constant  fCilisson  et  al  107S.  Moon  ef  at  1074)  as  well  as  electron  etieetiye 
mass  (Restorff  «7  al  1078).  In  these  cases  it  has  led  to  pood  agreement  with 
experimental  data,  l-'igure  1  shows  the  results  of  a  calculation  for  the  I' valley  electron 
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(ia,  ,ln»P,  yAs»  latticf-maichcil 
lo  InP  (.'ufvc  A  uses  binary 
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( Curves  R  aiul  C  use  values 
from  I  iiileiohn  ftal  ( W77a).  with 
(wo  ihffeteni  retHuteJ  values  for 
InP  V  ^  Kestorff  w<j/  (1^781 


effective  mass  in  Gai  ,In,l’i  ,ASv  as  a  function  of  As  composition  for  layers 
lattice-matched  to  Ini’.  The  three  different  curves  represent  calculations  based  on 
different  values  of  binary  material  parameters.  However,  for  y  r-O-.S  the  effective 
mass  does  not  strongly  depend  on  the  uncertainly  in  the  InP  effective  mass  (y  =0) 
fMaloncy  and  Frey  1077).  While  the  more  accepted  value  of  effective  mass  for  InP  is 
0  ()8  mo.  another  value  of  ()-0b7  m,,  has  been  used  in  previous  transport  calculations 
(Rode  1070)  and  is  shown  in  figure  1  for  completeness.  Also,  a  very  recent 
measurement  of  effective  mass  in  Ga,,  47lno  nAs  has  been  reported  as  ()-04 1  mo.  and 
the  interpolation  procedure  illustrated  in  figure  1  has  been  questioned  (Pearsall  fra/ 
1070).  This  measurement  was  made  at  near-helium  temperatures  from  Shubnikov- 
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dc  Haas  oscillatuin  ainpliliidos,  iii  the  saiiK'  maiiiK'i  as  Iho  iiicasiircmcnts  of  Rcstorif 
fi  al  (1‘>7S).  Ai  siii.h  loss  u-iups  i atiiics  the  coiiduelioii  band  density  of  states  is  well 
below  ltt"’em  and  one  woiikl  expeet  de)>enerate  statisties  to  be  prevalent  under 
the  tneasurement  eoiuhtunis  I'luis.  this  measured  elfeetive  masseould  be  larj^er  than 
the  band-edne  etfeetise  mass  shown  in  linure  I  ilue  to  non-parabolie  band  elleets.  In 
faet.  on  anothei  nioie  heasily  doped  (iad  uIiim  mAs  speeimen.  an  elfeetive  mass  ol 
1)1)74  ntii  was  tneasuted  This  mereasv  was  attributed  t»>  non-parabolie  band  etfeets 
(I’earsall  cf  a/  I  ‘)7'))  In  addition,  one  woulil  expeet  a  deerease  in  elfeetive  mass  with 
inereasm^  temperatuie.  whieh  wouhl  tend  to  reduee  the  helium  temperature  value  of 
t)  l)41  ni„  to  a  value  possibly  more  eonsistent  with  the  value  of  m,i  shown  in 

tifzure  1  for  Cla.i  47lni)  mAs.  The  expeiimental  elfeetive  mass  data  and  the  ealeulated 
values  shown  in  figure  1  appeal  to  be  eonsistent  at  this  point  in  time. 

In  general,  the  Monte  Carlo  method  is  dillieult  to  apply  with  good  statistical 
aceuraey  at  low  fields  and  low  temperatures.  However,  this  problem  ean  be  overcome 
to  a  large  extent  by  using  statistieal  estimators  for  the  dilfusion  coetiicient  (Canali  I’f 
a/  1975).  In  addition,  the  present  programme  has  been  further  modified  to  include 
the  effects  of  heavy  doping  as  well  as  to  use  crossed  electric  and  magnetic  fields  for 
the  computation  of  Hall  mobility,  since  it  is  Hall  mobility  and  not  drift  mobility  that  is 
usually  measured.  The  programme  details  will  be  presented  elsewhere.  Heavy 
doping  effects  can  be  important  in  Ga,  ,ln.P|  ,.Asv  alloys  having  small  effective 
masses  such  as  shown  in  figure  1.  For  GaAs,  with  m*  =  0-()67  mo,  the  conduction 
band  effective  density  of  states  is  about  4  x  lO'^  cm"’  at  300  K.  From  the  values  in 
figure  1,  a  number  of  1  -65  x  lo”  cm  '  for  the  density  of  states  of  a  quaternary  with 
y  =  0-7  is  calculated.  As  the  temperature  is  decreased  to  77  K,  the  conduction  band 
effective  density  of  states  is  reduced  to  5  X  lo"’cm  '  for  GaAs  and  to  2- 1  x  U)'*’ cm 
for  the  quaternary.  Thus  one  could  expect  degenerate  statistics  to  be  important  in 
determining  the  quaternarv  transport  properties  for  presently  achievable  doping 
levels. 

The  Monte  Carlo  melhrul  has  the  additional  feature  that  very  detailed  transport 
and  scattering  factors  can  be  included  in  the  calculations.  These  factors  include  such 
band  structure  effects  as  p-state  mixing  and  band  non-parabolicity  (Fawcett  et  al 
1970).  Also,  the  Monte  ('arlo  method  takes  into  account  thermal  stimulation  al 
temperatures  above  300  K  of  carriers  into  higher-lying  conduction  band  minima 
(Rode  1975).  If  care  is  taken  in  assuring  statistically  accurate  estimates  in  the  Monte 
Carlo  method,  then  the  technique  should  be  accurate  and  very  complete,  and  as 
convenient  as  techniques  based  on  classical  transport  analysis  and  the  use  of  other 
numerical  methods  (Rode  1975). 

4.  Discussion 

The  use  of  the  Monte  Carlo  method  for  calculating  Hall  mobility  was  first  examined 
by  considering  data  for  GaAs,  since  it  was  felt  that  values  for  all  important  transport 
coefficients  and  material  parameters  were  reasonably  well  known  (Rode  1975, 
Littlejohn  et  al  1977b.  Kratzer  and  Frey  1978).  Figure  2  shows  Hall  mobility  as  a 
function  of  free-electron  concentration  for  GaAs  at  300  K  for  a  compensation  ratio 
of  1  (that  is,  (N‘'  +  N  )/n  =  1.  where  N*  is  the  ionized  donor  density,  N  is  the 
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ionized  acceptor  density  and  n  is  the  tree-electron  density  ).  Shown  for  comparison  in 
this  figure  is  a  range  of  experimental  data,  taken  for  convenience  from  the  last  two 
proceedings  of  this  conference  and  the  standard  data  from  S/e  (1969).  Care  was 
taken  to  record  only  data  reported  as  Hall  mohility  data.  Also  shown  are  the  results 
of  a  calculation  for  Hall  mohility  based  on  an  iterative  integral  method  (Rode  1975). 
There  are  some  discrepancies  between  these  two  approaches,  as  can  be  seen  in  this 
figure.  Tor  low  carrier  concentrations  the  iliscrepancies  (about  1.5%)  can  be  attri¬ 
buted  to  slightly  ditferent  values  of  bulk  material  parameters  used  in  the  twi> 
calculation  methods.  The  mateiial  parameters  useil  here  in  the  Monte  Carlo 
calculation  have  been  chosen  baseil  on  both  low-licUl  mobility  values  and  on 
high-field  tr;itispr>rt  properties  (I  ittle|ohn  i7  nl  1977b,  Krat/er  and  Trey  1978). 
However,  the  Monte  Carlo  method  will  yield  Hall  mobilities  tit  low  carrier  concen¬ 
trations  ( s.  lO'^cm  ')  well  above  9()()()cm'  V  '  s  '  if  slight  adjustments  in  bulk 
parameters  are  made,  and  certainly  the  sensitivity  of  the  Hall  mobility  in  pure  GaAs 
to  slight  variations  in  material  parameters  is  well  recognized  (Rode  1975).  At  high 
carrier  concentrations  (>!()''’ cm  ')  the  deviations  between  these  two  compu¬ 
tational  methods  are  much  more  subtle.  They  are  due  to  the  way  in  which  the  ionized 
impurity  scattering  rate  is  corrected  for  degeneracy  in  the  Monte  Carlo  method,  and 
this  will  be  a  subject  for  a  future  paper.  However,  in  spite  of  the  differences,  the 
authors  believe  that  the  Monte  Carlo  method  is  an  accurate  technique  for  calculating 
Hall  mobility  which  gives  agreement  with  experimental  data  over  a  wide  range  of 
carrier  concentrations. 

Figure  .5  shows  experimental  Tfall  mobility  data  against  temperature  for  pure 
GaAs,  along  with  the  calculated  temperature  dependence  of  mobility  based  on  both 
the  iterative  integral  and  the  Monte  Carlo  method.  Here  the  Monte  Carlo  method 
gives  an  excellent  lit  over  the  temperature  range  1 50-600  K,  while  Rode  (1975)  has 
pointed  out  that  the  iterative  integral  lies  about  10%  above  the  experimental  data  in 
this  range  Whde  we  have  not  done  so  m  this  siuilv,  it  should  be  pointed  out  that  the 
Monte  C.irlo  method  can  be  exiended  to  icmpctaluics  below  10  K  (Canali  e’f  u! 
1975).  Also,  olhei  GaAs  malcii.il  gtown  in  oui  laboiaioiics  has  been  successfully 
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studied  by  the  Monte  Carlo  method,  and  the  general  use  of  the  method  appears  to  be 
well  justified. 

Shown  in  Figure  2  is  a  calculation  vtf  the  Hall  mobility  against  electron  concen¬ 
tration  at  300  K  for  Gap  :7lnn nPu  4AS11 ,,  based  on  the  Monte  Carlo  metht)d  and  the 
material  parameter  interpolation  procedure  described  previously  (Littlejohn  el  al 
iy77a).  At  the  present  lime,  to  the  authors'  knowledge,  the  largest  value  of 
room-temperature  Hall  mobililv  for  any  i|ualernary  is  abrnil  6000  em’  V  '  s  '  with 
an  electron  concentration  of  about  7x  lo''em  '  (Houston  el  ol  1078).  F'igure  4 
shows  a  summary  of  typical  Hall  mobility  agaitisl  temperature  for  the  i|ualernary 
(and  the  InP  lalliee-tnaleherl  ternary  (ia,,  .i;ln„  siAs)  taken  in  our  labrrralory  and 
reported  in  the  literature  (Hrnislon  el  til  107.S,  Takeda  el  til  1076).  In  general,  the 
room-temperature  Hall  mobilities  seem  to  lie  in  the  range  3000-4000  cm’  V  '  s  ', 


Figure  4.  Typical  Hall  mobility  data  against 
temperature  for  Gai  ,ln,Pi  yAs,.  The  curves 
show  data  from  Houston  et  at  (IV78I,  Takeda 
et  al  ( l*J7h),  and  from  this  study. 


244 


■V/  .4  l.illlcjohn.  R  A  Sthller.  /  //  (llisson  anti  J  R  Htiiiu’r 


with  clcctri>n  conccnlrjitioiis  of  (2-4)  x  lo"’  cm  Inlcicstiii>;lv  enough,  the  highest 
nu>bilitics  ;mil  lowest  clfctroii  coiKViilrations  coiisistc'iitlv  seem  to  occur  near  the 
ternary  hounilaries  lor  the  i|uaternarv.  that  is,  for  either  relatively  low  As  or  P 
compositions  This  is  sun^estive  of  the  inihienee  of  alloy  sealteiin>;  (Littlejohn  I'l  til 
l‘>7X). 

I'o  investiftate  the  >;eneial  piohlem  of  eai  i  iei  eompens.ition  aiul  alloy  seatterinp  in 
the  quaternary  (iai  .  In,  I’l  As,,  the  Monte  (  ai  lo  method  was  used  tostiuly  several 
typical  samples,  l  ifture  illustrates  the  situation  when  an  attempt  was  made  to  lit  the 
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Figure  5.  Hall  mobility  against  temperature 
for  (ia<i  \\no  7^0  46Asi)  m.  The  dotted  points 
arc  experimental  data  while  curves  A  and  R 
represent  Monte  ('arlo  calculations  with  no 
carrier  compensation  and  no  alloy  scattering, 
respectively 


data  using  either  carrier  compensation  with  no  alloy  scattering,  or  alloy  scattering 
with  no  carrier  compensation.  First,  the  alloy  scattering  potential  (Littlejohn  et  al 
1978)  was  varied  while  keeping  the  compensation  ratio  equal  to  l(i,e.  N*  +  N~  =  n) 
until  agreement  was  obtained  between  the  measured  and  calculated  Hall  mobility  at 
room  temperature.  Then  the  temperature  dependence  of  Hall  mobility  was  cal¬ 
culated,  resulting  in  curve  A  of  figure  .S  which  is  in  very  poor  agreement  with  the  data. 
Curve  B  results  by  removing  alloy  scattering  from  the  calculation.  First,  the  ionized 
impurity  density  was  varied  until  agreement  was  achieved  with  the  measured  Hall 
mobility  al  rinmi  temperature,  while  keeping  the  free-earrier  density  equal  to  its 
measured  riKim-temperature  value.  Curve  B  is  the  resulting  temperature  variation  of 
Hall  mobility  with  only  carrier  compensation  present.  Again  the  agreement  with 
experimental  values  at  other  temperatures  is  very  ptror,  and  it  was  concluded  that 
neither  alloy  scattering  nor  carrier  compensation  could  account  for  the  measured 
temperature  dependence  of  Hall  mobility. 
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Next,  an  attempt  was  maile  to  linj  a  eonsistent  set  of  values  for  holh  the  alloy 
scatterinj;  potential  anil  the  ioni/eil  impurity  density  whiili  would  together  pive  a 
good  fit  to  the  temperature  dependence  of  Hall  mobility.  Vo  do  this  two  experimental 
points  were  chosen,  one  at  .4(H)  K,  and  one  at  about  150  K.  A  detailed  parameter 
study  was  made  by  varying  the  alloy  scattering  potential  and  ioni/ed  impurity  density 
until  a  set  of  values  for  these  twi>  parameters  whieh  agreed  with  the  measured  values 
of  both  Hall  mobility  and  electron  concentration  at  the  two  temperature  extremes. 
Then  the  complete  Hall  mobility  against  temperature  curve  was  calculated  using 
these  values.  The  results  are  shown  in  ligure  (i  I'hc  Monte  Carlo  estimates  lie  within 


Figure  6.  Experimental  data  and  Monte 
Carlo  calculations  of  Hall  mobility  against 
temperature  for  Ga,,  ,lno  7P0  4sAs,i  »4 
showing  titted  material  parameters.  .V,,  = 
77  X  to'- cm  ’,.VA  =  5-5xl0'''cm  '.Ea  = 
0  7h  cV:  compensation  ratio  =  3  4. 


the  broken  curves  (several  estimates  were  maile  at  each  temperature)  while  the  bars 
indicate  a  5%  variation  in  the  experimental  data  at  the  two  temperature  extremes.  It 
should  be  pointed  out  that  the  fits  are  facilitated  by  the  fact  that  the  electron 
concentration  does  not  vary  significantly  for  these  quaternary  samples  between 
100-300  K..  The  compensation  ratio  of  3  ‘)  is  very  close  to  a  sample  with  a  very 
similar  temperature  dependence  of  mobility  studied  by  Houston  el  al  (1^)78). 

The  calculations  arc  more  sensitive  to  the  ionised  impurity  density  than  they  are  to 
the  alloy  scattering  potential.  I'hiis,  it  is  not  feasible  to  draw  many  conclusions  from 
the  actual  value  of  the  alloy  scattering  parameter  shown  in  figure  h.  The  value  is 
larger  than  the  proposed  theoretical  values  for  three  different  models  (Little  john  eltd 
1978)  by  as  much  as  20  to  100%.  There  is  another  scattering  mechanism  which  has 
the  same  temperature  dependence  on  mobility  as  alloy  scattering,  which  is  often 
called  space-charge  scattering  (Weisberg  1962).  If  this  mechanism  were  also  present 
in  these  samples  it  would  be  difficult  to  detect  from  Hall  measurements  alone,  but 
could  effect  the  determination  of  the  alloy  scattering  pcitential  by  the  technique  used 
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seeond-phase  ilefeels 

Other  samples  with  similar  Hall  mohility  temperature  ilepeiulenees  have  been 
studied  with  very  similar  results  as  niven  in  tiRure  h  In  all  eases,  the  alloy  seatterin>; 
potential  required  to  lit  the  data  was  laiper  than  the  theoretieal  value.  Also,  (or  a 
given  compensation  ratio,  a  weaker  temperature  dependence  of  mohility  is  indica¬ 
tive  of  a  larger  effect  due  to  alloy  scatter  ing.  This  is  not  inconsistent  with  figure  5. 
since  there  the  room-temperature  mobility  was  used  as  the  target  value  to  be  fitted  by 
varying  only  a  single  parameter;  either  the  compensation  ratio  or  the  alloy  scattering 
potential. 

It  is  interesting  to  consider  those  samples  w  ith  a  stronger  temperature  dependence 
in  this  temperature  range,  ^■igurc  7  shows  the  results  for  a  (ia„  ailuo  mAs  sample. 


T- 

i. 


C 


Tomperoture  IK) 


Figure  7.  Experimental  data  and  Monte 
C'arlo  calculations  of  Hall  mobility  against 
temperature  for  Ga,,  47lno  <  »As  showing 
fitted  material  parameters.  .Vn=»l5x 
tO'^cm  iVA=  lOx  I0''’cm  = 

0-4  eV';  compensation  ratio  =  4-8. 


which  is  a  ternary  in  this  quaternary  system  having  a  Hall  mobility  at  room 
temperature  larger  than  8000  ctn  V  '  s  '  (Morkoc  1478.  Takeda  1478).  I'or  this 
ease,  the  temperature  dependence  of  Hall  mobility  could  be  lit  by  the  Monte  Carlo 
method,  but  with  an  alloy  scattering  potential  essentially  the  same  as  that  predicted 
by  the  theoretical  model.  The  influence  of  alloy  scattering  in  this  ternary  is  less 
important  in  determining  the  Hall  mobility  than  it  is  in  the  quaternary.  1  his  fact  is 
somewhat  surprising  based  on  previous  results  for  the  scattering  potentials  in  both 
ternary  and  quaternary  alloys  (l.ittlejohn  el  al  1478). 


( 'timer  t  t>ini>t'ii\tilion  iiiul  tillin’  yi  tilU’riiifi  in  (hi ,  JiiJ’i  ,.4  v, 
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5.  Conclusions 

Alloy  scattermk:  appears  ir'  play  an  iiiiporiaiit  role  m  the  low-tiekl  transport 
properties  of  Gai  ,  In,  I’l  ,  As,.  This  eonelusion  is  based  i>n  the  data  and  ealeulations 
presented  in  this  paper,  as  well  as  on  other  published  Hall-etleet  data.  However,  the 
accuracy  of  the  ealeulations  is  (uedieated  on  the  material  paiameter  interpolation 
procedure,  which  is  presently  renuireil  for  anv  transport  ealeulations  in  this  quater¬ 
nary.  Mistorieallv.  the  Monte  (  arlo  method  has  become  more  useful  in  studying 
transport  properties  m  eompounil  semiei'iuluetors  as  accurate  material  parameters 
become  available  Hie  iletei inmatioii  of  the  basic  material  parameters  should 
be  an  important  part  of  progi  aninies  which  attempt  to  obtain  high-purity 
(5ai  ,  In, I’l  ,.  As,  quaternary  alloys 
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A  nethod  of  including  the  ^^auli  Exclusion  Principle  in  the  Monte 
Carlo  procedure  is  presented.  The  procedure  is  then  used  to 
calculate  the  Hall  nobility  and  velocity- field  characteristics 
tor  GaAs  with  impurity  concentrations  up  to  10>»  cm-3. 


Introduction 

The  Monte  Carlo  technique,  first  described  by  Kurosawa'  1  ]  and 
then  generalized  by  Fawcett  ct  al' 2  ],  has  become  a  very  important 
tool  in  predicting  semiconductor  material  characteristics.  while 
it  has  been  used  to  analyze  many  different  materials,  it  is 
usually  restricted  to  nondegenerate  cases.  To  remove  this 
restriction  from  the  usual  Monte  Carlo  procedure  electron- 
electron  scattering  and  the  Pauli  Exclusion  Principle  must  be 
included.  Some  attempts  to  include  the  former  have  been  reported 
[1-6]  and  one  to  include  the  latter  has  been  reported[7]. 

In  this  paper,  one  method  of  incorporating  the  Pauli  Exclusion 
'’rinciple  in  the  .Monte  Carlo  procedure  is  presented  which  is 
valid  at  low  fields.  The  method  presented  here  differs  from  that 
described  by  S.  Posi  and  C-  .1acoboni[7]  in  that  (1)  an  estimate 
of  the  probability  of  occupancy  is  used  instead  of  accumulating 
this  information  via  a  k-snacc  mesh  during  the  simulation  and  (2) 
ionized  impurity  interaction  is  treated  as  a  special  case.  This 
method  is  used  to  calculate  the  low  field  Hall  mobility  of  GaAs. 
Ar  ex^rension  of  this  method  is  then  used  to  deteraine  the 
velocity-t iel d  characteristics  of  GaAs.  Both  calculations  are 
done  tor  both  nondegenerate  and  degenerate  cases. 

The  Monte  Carlo  model  used  in  this  work  is  the  same  as  that 
described  elsawhorr  [8,9]  except  the  present  model  includes  p- 
G^ate  mixing.  Briefly,  it  inconoratos  a  three  valley  model  with 
nonpa ra bolici t y ,  and  acoustic  phonon,  optical  phonon,  ioniztd 
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impurity,  polar  optical,  piczofilcctr  ic,  and  intervalley 
scattering  processes.  The  raa^^erial  parameters  used  for  GaAs  are 
the  sane  as  those  presented  in  [  10  ]  and  are  repeated  in  Table  1 
tor  completeness. 


Corrections  to  Scattering  Rates 

The  Pauli  Exclusion  Principle  must  bo  considared  in  xodelling 
degenerate  semiconductors.  To  include  the  Pauli  Exclusion 
Principle  in  the  Monte  Carlo  procedure  described  in  the  previous 
section,  we  multiply  the  scattering  rate  Sj^(k,k')  between  two 
wave-vector  states  k  and  k*  by  the  factor  (l-f(k')).  Here  f(k') 
IS  the  probability  the  final  state  k'  is  occupied  so  that 
(1-f(k*))  is  the  probab  lity  the  state  is  unoccupied.  Thus,  the 
total  scattering  rate  from  the  state  k  due  to  the  n'th  process  is 

X  (k)  =  /dk'  S  (k,k')(l-f  .  (1) 

n  n 


For  low  fields,  the  occupation  probability  distrubition  f (k) 
is  near  equilibrium,  so  that 

f(lt)  =  f^(t)  (2) 

where  f  ^  (k)  is  the  equilibrium  occupation  probability 
distribution  which  is  the  Fermi-Dirac  distribution 


f  (k) 


(3) 


In  the  work  reported  here,  the  relationship  between  the  energy  E 
and  wave-vector  k  ir.  the  i*th  valley  (i=1,2,3)  is  taken  to  be  [2] 


^  .  ,,(E) 


where  a  is  the  nonoarabolicit y  and  m*  is  the  effective  mass. 


To  correct  the  scattering  rates  as  in  Eg.  (1),  only  the  Fermi 
energy  need  be  calculated.  This  is  done  by  solving  (numerically) 
the  equation  [11] 


V  u  *3/2 
V  Z  M  ,  m* 

2.  3  .  ,  ci  i 

71  h  i=l 


A^(E) 


dy^CE) 


f  (E)dE 
o 


where  Md  is  the  number  of  equivalent  minima  in  the  i’th  valley 
and  n  is  the  density  of  conduction  band  electrons. 


Pigurp  1  shows  thp  scattorinrj  rate  for  polar  optical  emission 
with  and  without  the  correction  factor  in  the  valley  for  an 
irpirity  concentration  of  lO***  cm”’  and  a  temperature  of  300  K. 
In  general,  the  correction  factor  effectively  prohibits 
scatterings  that  would  result  in  a  final  energy  less  than  the 
Fermi  energy. 

The  correction  factor  is  included  in  all  scattering  rates 
except  ionizel  impurity.  in  ionized  impurity  creates  a  localized 
perturbation  on  the  electric  field.  The  wave-vector  k  of  an 
electron  influenced  by  this  pt.'rturbing  field  is  not  changed 
abruptly  to  k*  but  is  gradually  changed  over  the  time  the  carrier 
is  influenced  by  the  perturbing  field.  Thus,  even  though  an 
ionized  impurity  interaction  is  treated  in  the  nonte  Carlo 
procedure  as  an  instantaneous  scattering  from  wave-vector  k  to 
k’,  the  correction  factor  (l-t(k’))  seems  inappropriate. 


Deteraination  of  Low  Field  Hall  Kobility 

To  calculate  the  Hail  mobility  using  the  Monte  Carlo 
procedure,  it  is  necessary  to  include  a  magnetic  field.  For  this 
case  the  time  derivative  of  the  wave-vector  k  is 


3  t 


^  (c  +  V  X  B) 

n  n 


(6) 


where  c  is  the  electric  field  vector,  B  is  the  magnetic  field 
vector,  and  w  is  the  carrier  velocity. 

flere  he  electric  field  is  taken  to  be  in  the  z  direction  and 
the  magnetic  field  in  the  x  direction.  Vith  this  configuration 
the  low  field  Hall  mobility  is  shown  in  the  appendix  to  be  given 
by 


y  z 


where  D  y,  and  are  the  diffusion  coefficients  in  the  x,  y, 
and  z  directions,  respectively,  and  P  is  the  magnitude  of  the 
magnetic  field.  Since  P  and  n  are  equivalent  for  low  fields, 
the  estimate  of  u  may  be^ improved  by  using 


yz 


(8) 


where 


3 


(9) 


D 

yz 


=  i  (D  +  D  ) 


Using  "gs.  (3)  and  (3)  in  tha  Monte  Carlo  progra>  with  a 
magnetic  tieU,  the  Hall  mobility  versus  impurity  concentration 
was  calculated  for  CaAs  with  corrections  to  the  scattering  rates 
as  described  earlier  (m  correction  to  the  rate  for  ionized 
impurity  scattering).  Figure  2  shows  a  simple  flow  chart  of  the 
program  used.  The  calculation  was  also  carried  out  for  two  other 
cases:  (1)  no  correction  to  the  scattering  rates,  and  (2) 
correction  to  all  the  scattering  rates  (including  ionized 
impurity).  The  results  are  given  in  Pig.  1.  Also  shown  is  the 
e xp<' r  1  men t a  1  data  obtained  by  Szn[12].  Sood  agreement  with  ^ze's 
experimontil  lata  is  obtained  by  using  the  correction  to  all  the 
scattering  ratt‘s  except  ionized  impurity,  while  the  other  two 
cases  give  rather  poor  agreement,  supporting  the  special 
treatment  of  the  ionized  impurity  interaction. 


Calculation  of  Velocity-Field  Characteristics 

The  method  used  above  to  account  for  the  Pauli  Exclusion 
Principle  in  the  calcula*-ion  of  low  field  Hall  mobility  must  be 
modified  for  the  calculation  of  velocity- field  characteristics 
tor  degenerate  semiconductors.  The  method  adopted  here  is  to 
assume  the  probability  of  occupancy  in  the  i'th  valley  {i=1,2,3) 

IS 


f^(E)  = 


1+c 


(10) 


Thus  a  quasi  Fermi  level  and  quasi  electron  temperature  Tei 
art*  associated  with  the  carriers  in  the  i'th  valley.  The 
correction  to  the  scattering  rates  in  the  i'th  valley  is  then 
based  on  the  probability  of  occupancy  given  in  Eg.  (10). 


Incorporating  tnese  new  features  in  the  Tonte  Carlo  procedure 
requires  an  iterative  solution  for  the  carrier  temperatures  and 
the  Fermi  erergies-  Initial  values  are  determined  by  assuming 
that  (1)  the  carrier  temperatures  are  all  equal  to  the  lattice 
temperature  and  (2)  the  Fermi  energies  are  all  equal.  Then,  Eg. 
(S)  can  be  solved  tor  the  Fermi  energy.  Using  these  initial 
values,  the  scattering  rates  are  corrected  and  the  carrier  is 
tracked  tor  some  specified  number  of  interactions.  The  density 
of  carriers  and  the  average  energy  S.  in  each  valley  are  then 
calculated.  Then  ^ 


n , 
1 


/2 


IT  h 


M  .  m  .  * 
cl  i 


3/2 


'TOPI 


A^(E) 


dY.(E) 


Cl 
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iin  J 
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M  ,  m  * 
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2,  3 
n  h 


3/2 
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roi'i 

E»  >  j  (E) 


dYjCK) 

ciE 


lj(E)  dE 
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aro  usod  to  obtain  now  o:;t  i  na  to*;  of  tho  quasi -o  1  oct  i  on 
t  ora  jx*  ra  t  uros  and  tho  quasi-Formi  onorqios.  Tho  now  cairior 
t  om  po  ra  t  ur  o  s  and  Forrai  oiu*iqios  aro  t  hi'n  usod  to  obtain  now 

corroct ion  factors  {1-t  (k'))  and  now  scattorinq  rates  and  tho 

process  is  repeated.  This  proc(*ss  of  corroctinq  scatterinq 

rates,  tiarkinq  the  earlier,  and  calculating  new  ostinates  of 
carrier  temperatures  and  Fermi  onorqios  is  continued  until  a 
so  1 1 -con s 1 st ent  solution  is  reached,  i-o-,  until  succossivo 

iterations  give  little  change  in  Epj  and  .  Figure  U-  shows  a 
tlow  chart  tor  the  iterative  procedure  lor  a  single  valley. 


Using  tho  method  outlined  above  in  the  Monte  Carlo  procr'duro, 
the  VO  loci  t  y- f  i  e  1  i  ch  a  r  ac  t  or  i  st  i  cs  of  r^aAs  were  calculated  for 
I'lpurity  concentrations  of  10‘*,  10*",  and  fO*"*  cb“>.  Tho 

results  are  given  in  Fig.  S.  For  comparison  tho  velocity-field 
characteristics  for  GaAs  calculated  from  tho  Monte  Carlo 
procedure  without  incluJing  tho  Pauli  Exclusion  Principle  are 
shown.  For  the*  nondogenera  t  e  case,  impurity  concentration  of 
10“*  cm~^,  Yhe  two  curves  are  in  good  agieoment,  while  for  the 
two  legenerate  caj-.es  the  difference  becomes  more  pronounced. 
This  result  is  expected  since  the  Pauli  Exclusion  Principle  is 
important  only  tor  ttu'  decu*n»'ratc  cases. 


A  pppnd I X 


An  «‘xpr!>ssion  ti)r  thr  Il-ill  mobility  in  t«'rm 
coot  t  iciont  s  c<in  l)o  ot’tainod  t  rom  tho  cut  rent 
tor  low  tii*13s 


of  t  ho  lift usion 
d  (Ml  si  t  y  p(]  ua  t  i  on 


in  terms  of  which  Fq.  (A1)  can  be  written 


IS  of  the  form 


he  solution  to  Fq 


Substitution  of  Eq 


where  Cj,  C  2,  and  C3  are  to  be  determined 
(A5)  into  the  right  side  of  F.q.  (A4)  yield 


Fquiting  coefficients  between  Fqs.  (A5)  and  (A6)  gives 


Fq.  (Ah)  becomes 


Vn)]B  +  W,,(c-t— 


[  B*(.  (4— 


where 


BL-C. 


For  the  car.e  where  the  raaejnotic  field  is  in  the  x  direction 
and  the  electric  field  is  in  the  z  direction 

Vn-B=g^B.  (A  9) 


cxB  =  -  tBp, 


(A10) 


II  n  nSn- 

VnxB  =  B  a—  y  -  B  37 


(All) 


and  Fq.  (A7)  becomes 


-t  3n„  ^  n _  /3  n  n  ^  " 

57  «  “72  ^7  57»' 


qg  n(u  B) 
n  H 


2  2  c 


n 

y  + - E  z 


y  -TT  (A12) 

1+Uh  b^  b 

.  /in  _  D  ?Ji\  " 

2  2  S  z  3y^  • 

1+uu  B  ■' 

NOW,  assuminq  the  gereratioi.  and  recombination  rates  are  equal/ 
the  continuity  equation  for  electrons  is 


is  .  i  V  .  3 
3  t  q  n 


(All) 


Substitutinn  Fq.  (M2)  i.uo  Fg.  (A13)  yields 


Lq  =  ,  ^  ^  -JJl  ^  U  ^  ^ 

3t  z3z  y3y  X3^2  >.^^2 


(A1«) 


where 


z  I.  2„2  ’ 

1+Uh  b 


(A15) 


-U  E(g  B) 
n  H 

''v  "  - ^ 


(A16) 


(^17) 


D  =  I) 
X  n 


D  =  D 

y  z 


D 

n 


2  2 

IH-U,,  B 


(A18) 


Pinaily,  combining  Eqs.  (A  17)  and  (A18)  yields 


which  is  the  desired  result. 


(A19) 
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Tilblo  1.  CtdAs  r.atprial  parameters  us'*d  in  calculations. 


"ulk  aaterial  parameters 
Parameter 

Value 

Lattice  constant 

5.642X10 

- »  cm 

Density 

5.36  g/cm^ 

Electron  affinity 

4.07  eV 

Piezoelectric  constant 

0.16  C/m 

? 

LO  phonon  energy 

0.03536 

eV 

Sound  velocity 

5. 24X  10* 

o/sec 

Optical  dielectric  constant 

10.92 

Static  dielectric  constant 

12.90 

Valley-dependent  material  parameters 

Conduction-band  valley 

Parameter 

r(OOO) 

X (100) 

L(111) 

Acoustic  lef.  pot.  (eVJ 

7.0 

9.27 

9.  2 

Effective  mass  (m  /m  ) 

0.063 

0.  58 

0.222 

Nonparabolicit  y  (eV-») 

0.610 

0.204 

3.461 

Energy  band  gap  (eV) 

1. 43P 

1.961 

1.769 

(relative  to  valence  band) 
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INFLUENCE  OF  CENTRAL  VALLEY  EFFECTIVE  MASS  AND 

ALLOY  SCATTERING  ON  TRANSIENT  DRIFT  VELOCITY  IN  Ga,  In  P  As  * 

1-x  X  1-y  y 

M.  A.  Littlejohn,  L.  A.  Arledge,  T.  H.  Glisson,  and  J.  R.  Hauser 
Electrical  Engineering  Department 
North  Carolina  State  University 
Raleigh,  NC  27650 


ABSTRACT 


The  instantaneous  ensemble-averaged  and  time-averaged  velocity  of 

electrons  are  presented  for  Ga  T^In  -,-P  ,As  ,,  and  the  results  are 

.2.1  , 1  i  .  ^  .0 

compared  to  those  for  GaAs.  The  results  indicate  that  both  the  effective 
mass  in  the  F  conduction  band  and  alloy  scattering  are  very  critical  in 
determining  the  transient  (or  velocity-overshoot)  characteristic  as  well 
as  the  static  velocity-field  characteristic  for  this  quaternary  material. 


*This  work  was  supported  by  the  Office  of  Naval  Research,  Arlington, 
Virginia,  U.S.A. 


In  a  previous  publication  [1]  the  static  (time-averaged)  velocity 


field  characteristic  of  the  Ca,  In  1’ 

1-x  X  1 


presented.  This  cnaracter  ist  ic  wa.s  calculated  using  a  monte  carlo 


simulation  of  the  electron  transport.  The  quaternary  material  properties 


used  in  the  simulation  were  obtained  by  an  interpolation  procedure 


based  on  experimental  and  theoretical  material  parameters  for  the  binary 


constituents  of  the  quaternary.  This  interpolation  procedure  for  a 


given  material  parameter  results  in  a  continuous,  smoothly-varying 


function  in  the  (x,y)  compositional  plane,  which  reduces  to  the  correct 


compositional  dependence  along  the  ternary  boundaries  and  to  the  binary 


values  at  the  corners  of  the  compositional  plane 


There  have  been  very  few  experimental  investigations  of  the  basic 


material  properties  of  Ga  In  P.  As_,  and  thus  very  few  means  for 


confirmation  of  the  above-mentioned  interpolation  procedure.  Energy 


band  gap  measurements  by  photoluminescence  [2-4]  do  support  the  interpo 


lation  procedure.  However,  the  drift  velocity  is  affected  onlv  very 


slightly  by  the  value  of  the  band  gap.  Recently,  several  techniques 


have  been  used  to  measure  the  band  edge  effective  mass  [5-7],  and  this 


parameter  fundamentally  affects  the  transport  properties  of  any  semi 


conductor.  Some  of  these  measurements  do  support  the  interpolation 


technique  [5,6].  However,  other  experimental  data,  combined  with 


theoretical  calculations,  [7]  indicate  that  the  effective  mass  has 


substantially  less  "bowing"  than  is  predicted  by  the  interpolation 


formula  [8]  for  quaternary  compositions  which  are  lattice-matched  to  InP 


Very  recently-presented  data  ter  botli  electrnn  mobility  (9)  and  electron 

ettectlve  mass  (10)  turtlier  support  this  contention. 

The  purpose  ot  this  paper  is  to  examine  the  Inlluence  ol  central  (!') 

valley  etlectlve  mass  on  the  calculations  4)1  both  the  time-averaged  static 

dritt  vel4)cltv  vet  sus  electilc  tlcKl  Intensity,  itnd  the  ensemble-averajted 

ti.ittslent  veUx-.ltv  (veh'cltv  ov4'rshoot  )  c hariicter  1  st  1  c  lor  (la,  In  1’.  As  . 

1-x  X  1-y  y 

Also,  the  Impi'rtance  ol  alK'y  Si-atterlnp,  will  be  lurther  considered  (1). 

The  calculations  itre  presented  tor  one  composition  (x“0.7J,  y“0.t))  whicli 

is  latt ice-nuitched  to  Ini’  anil  has  a  room-temperature  band  gap  ol  0.96  eV, 

since  this  particular  compositi4)n  has  been  suggested  as  being  suitable 

lor  microwave  device  applications  [1]. 

Transient  velocity  calculations  h.tve  not  been  reported  for  the 

quaternary,  and  tlie  concept  of  velocity  overshoot  could  become  very 

imp4)rt.int  tor  devices  with  submicron  active  dimensions.  It  will  be 

slu)wn  In  this  paper  that  the  T-valley  etlectlve  muss  Is  of  basic 

imp4)rtance  in  determining  the  niiture  of  the  veloclty-4)vershool  characteristic 

ot  Oa,  In  r.  As  .  One  pr Imarv  re.ison  Is  that  this  quaternary  Is  a  material 
1-x  X  1-y  y 

which,  to  a  liirge  extent.  Is  di)minated  by  scattering  in  the  central  valley 
due  to  the  predictevl  large  I’-l.  Intervalley  energy  separation  (11). 

The  static  velocity-field  calculations  have  been  described  previously 
[1).  In  the  transient  analysis  procedure,  electrons  are  rjindomly  selected 
from  a  Mav-_-£;iii£.r|  velocity  distribution  with  T-300“K.  The  electric 


field  intensity  is  independent  of  position.  Each  electron  undergoes 


3 

scattering  as  it  would  in  a  single-electron  monte  carlo  procedure. 

During  the  drift,  the  instantaneous  velocity  and  position  of  the  electron 
are  recorded  at  specific  times  over  a  time  interval  generally 
long  enough  to  allow  the  particle  to  reach  a  steady-state  velocity.  This 
procedure  is  then  repeated  N  times  (N“5000),  and  the  instantaneous  ensemble 
averaged  velocity  and  distance  along  the  field  direction  are  calculated  as 
a  function  of  the  time.  The  data  presented  liere  will  be  in  the  form  of 
velocity  versus  distance  since  this  gives  a  more  qualitative  assessment 
of  the  persistence  of  velocity  overshoot  in  terms  of  physical  device 
dimensions . 

Figure  1  shows  the  static  velocity-field  characteristic  for 

Ga  27^*^  73^  4'^®  5  values  of  F-valley  effective  mass,  along  with 

17  -3 

a  curve  for  GaAs.  The  carrier  density  is  1x10  cm  .  The  effective 

mass  of  0.036  m^  for  the  quaternary  is  obtained  from  the  interpolation 

procedure  [1,8],  while  the  value  of  0.05  m^  is  chosen  based  on  the  more 

recently-measured  values  [7,10].  The  calculation  includes  alloy  scattering 

for  the  quaternary  and  assumes  a  random  alloy  with  no  sublattice  ordering 

[1,12].  The  increased  F-valley  effective  mass  decreases  the  drift 

2  2 

mobility  of  the  quaternary  from  7000  cm  /volt  sec  to  4900  cm  /volt  sec, 
which  is  very  close  to  the  GaAs  value.  The  peak  velocity  is  still  larger 
than  for  GaAs  at  this  doping  level,  although  the  threshold  field  is  also 
larger.  The  peak-to-valley  ratio,  with  the  valley  velocity  taken  at 
20  kV/cm,  is  nearly  the  same  for  GaAs  and  the  quaternary  with  m*(F)«0.05  m^ 
Also  shown  in  Figure  1  is  the  velocity-field  characteristic  for  the 
quaternary  with  m*(r)“0.05  m^  without  alloy  scattering.  Alloy  scattering 
is  more  detrimental  to  the  static  velocity-field  curve  for  this  value  of 


4 


effective  mass  than  for  the  value  of  0.036  m  [11.  Some  of  the  recent 

o 

experimental  transport  property  data  does  Indicate  that  a  scattering 
mechanism  having  a  temperature  dependence  of  mobility  the  same  as  alloy 
scattering  is  greatly  affecting  the  electron  dynamics.  The  problem  of 
alloy  scattering  in  the  quaternary  material  remains  as  a  major  barrier 
to  the  realization  of  many  of  the  potentially  useful  material  properties. 

Figure  2  shows  the  velocity  versus  distance  characteristic  for  the 

quaternary  both  with  and  without  alloy  scattering  for  effective  mass 

values  of  0.036  m  and  0.05  m  and  a  constant  field  of  E=10  kV/cm. 
o  o 

Shown  for  comparison  is  a  curve  for  GaAs,  which  is  in  very  good  agree¬ 
ment  with  results  calculated  by  other  authors  [13-15] .  Figure  3  shows 
similar  calculations  for  a  field  of  25  kV/cm.  There  are  several  points 
to  be  made.  For  an  effective  mass  of  0.036  m^  the  velocity-overshoot  is 
very  dramatic,  and  even  including  the  effects  of  alloy  scattering  it 
persists  over  distances  of  greater  than  approximately  0.35  pm.  For  an 
effective  mass  of  0.05  m^  the  velocity  overshoot  is  much  reduced  and  for 
small  distances  (0.3  pm)  the  Instantaneous  velocity  for  GaAs  is  greater 
than  that  for  the  quaternary  with  alloy  scattering.  It  should  be  pointed 
out  that  due  to  the  shift  of  the  threshold  field  in  the  quaternary  with 
m*(r)  =  0.05  m^  to  about  8  kV/cm,  one  would  expect  a  very  small  amount  of 
velocity  overshoot  for  fields  only  slightly  above  the  threshold  field 
[13-15],  as  is  the  case  in  Figure  2.  The  effect  of  alloy  scattering  on 
the  velocity  overshoot  for  m*(r)  »  0.05  m^  is  considerably  larger  than 
that  for  m*(r)  “  0.036  m^.  It  can  be  seen  that  both  alloy  scattering 
and  central  valley  effective  mass  are  predicted  to  have  a  significant 
effect  on  the  velocity  overshoot  characteristic  for  this  quaternary  alloy. 


5 


and  Lhe  actual  values  of  parameters  used  in  the  simulation  will  determine 

whether  or  not  Ga,  In  P,  As  is  far  superior  to  GaAs  or  is  only 
1-x  X  1-y  y  ^ 

comparable  to  GaAs. 

Recently,  it  has  become  apparent  that  new  physical  concepts  could 

be  involved  in  the  electron  dynamics  in  sub-micron  device  structures  [16] . 

Also,  it  is  obvious  that  for  real  devices,  such  as  in  the  channel  region 

of  a  micron-geometry  MESFET,  the  electric  field  intensity  is  not  constant 

and  varies  drastically  with  position  in  the  channel.  The  spatial  variation 

of  electric  field  will  play  a  significant  role  in  the  influence  of 

electron  dynamics  on  device  performance.  Thus,  calculations  such  as  those 

presented  here  and  elsewhere  [13-16]  can  only  be  taken  as  the  first  step 

in  assessing  the  magnitude  of  the  velocity-overshoot  effect.  By  comparing 

two  materials  in  the  manner  presented  here,  one  can  assess  the  relative 

performance  of  the  quaternary  in  comparison  to  GaAs.  Based  upon  the  work 

presented  here  there  still  appear  to  be  reasons  for  being  optimistic  about 

the  potential  of  Ga,  In  P,  As  as  a  microwave  material.  However,  as 
1-x  X  1-y  y 

these  results  demonstrate  any  definite  conclusions  must  await  further 
experimental  evaluation  of  the  exact  material  parameters  of  the  quaternary. 
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Figure  1.  Static  drift  velocity  versus  electric  field  intensity  for 
T=300K  and  Np=lxl0^^cin“^ .  The  notation  is: 

-  Ga  „_ln  ^.P  .As  ,  C-  GaAs 

A  -  in*(r)  =  0.036  with  alloy  scattering 

B  -  m*(r)  =  0.05  with  alloy  scattering 

B'-  m*(r)  =  0.05  without  alloy  scattering 


Figure  2.  Transient  drift  velocity  versus  distance  for  an  electric 
field  of  10  kV/cm  with  T=300K  and  Np=lxl0^^cm“^ .  The 
notation  is: 


A,A',B,B'  -  Ga  27^^  73^  4^®  5 


C-  GaAs 


A  -  in*(r)  =  0.036  with  alloy  scattering 

A'-  n)*(r)  =  0.036  without  alloy  scattering 


B  -  m*(r)  =  0.05  m 

c 

B'-  m*(r)  =  0.05  m 


with  alloy  scattering 
without  alloy  scattering 


Figure  3.  Transient  drift  velocity  versus  distance  for  an  electric 

field  of  25  kV/cm  with  T=>300K  and  Np=lxl0^^CDi“^ .  The  notation 
is  identical  to  that  of  Figure  2. 
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